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Abstract. 

The supersymmetric method is a powerful method for the nonperturbative evaluation of quenched 
averages in disordered systems. Among others, this method has been applied to the statistical 
theory of S-matrix fluctuations, the theory of universal conductance fluctuations and the microscopic 
spectral density of the QCD Dirac operator. 

We start this series of lectures with a general review of Random Matrix Theory and the statistical 
theory of spectra. An elementary introduction of the supersymmetric method in Random Matrix 
Theory is given in the second and third lecture. We will show that a Random Matrix Theory can 
be rewritten as an integral over a supermanifold. This integral will be worked out in detail for the 
Gaussian Unitary Ensemble that describes level correlations in systems with broken time reversal 
invariance. We especially emphasize the role of symmetries. 

As a second example of the application of the supersymemtric method we discuss the calculation 
of the microscopic spectral density of the QCD Dirac operator This is the eigenvalue density near 
zero on the scale of the average level spacing which is known to be given by chiral Random Matrix 
Theory. Also in this case we use symmetry considerations to rewrite the generating function for the 
resolvent as an integral over a supermanifold. 

The main topic of the second last lecture is the recent developments on the relation between the 
supersymmetric partition function and integrable hierarchies (in our case the Toda lattice hierarchy). 
We will show that this relation is an efficient way to calculate superintegrals. Several examples 
that were given in previous lectures will be worked out by means of this new method. Finally, 
we will discuss the quenched QCD Dirac spectrum at nonzero chemical potential. Because of the 
nonhermiticity of the Dirac operator the usual supersymmetric method has not been successful in 
this case. However, we will show that the supersymmetric partition function can be evaluated by 
means of the replica limit of the Toda lattice equation. 
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1. INTRODUCTION 



Random Matrix Theories are ensembles of matrices with random numbers as matrix 
elements. They were first introduced in physics in 1955 by Wigner [1] to describe the 
spacing distribution of nuclear levels. The mathematical problem of Random Matrix 
Theory is to derive eigenvalue correlations in the limit of large matrices. This problem 
was solved by Dyson [2] and Mehta and Gaudin [3, 4] who introduced orthogonal 
polynomials to perform the matrix integrals. About the same time the Anderson model 
[5] was introduced to describe the localization of wave functions in one-dimensional 
systems which led to the theory of disordered systems. Random Matrix Theory, with 
its main applications to nuclear physics, and the theory of disordered systems, with 
its main applications to solid state physics, developed independently until the mid 
eighties when it was realized that Random Matrix Theory is the zero dimensional 
limit of a disordered system [7, 6, 8, 9]. This was made possible in particular by the 
introduction of a field theoretical formulation of the theory of the disordered systems 
[10], and by the introduction of the supersymmetric method which made it possible 
to perform nonperturbative calculations. After these developments a rapid growth of 
Random Matrix Theory followed in the late eighties after the discovery of universal 
conductance fluctuations [11] and applications of Random Matrix Theory to quantum 
gravity [12, 13]. Both applications greatly enhanced the intellectual breath of the subject. 

Random Matrix Theory has many applications. We already mentioned the correlations 
of nuclear levels [14]. It was realized in 1985 [15] that the necessary ingredient for levels 
to be correlated according to Random Matrix Theory is not the complexity of the system 
but rather that the system is classically chaotic. This led to many applications in the field 
of Quantum Chaos. A second application of Random Matrix Theory in nuclear physics 
is the description of 5-matrix fluctuations . This problem was solved by means of the 
supersymmetric method [9] and the same approach was later applied to the theory of 
universal conductance fluctuations [16]. 

A new class of Random Matrix Ensembles was introduced in the early nineties 
[17, 18] to describe the eigenvalue correlation of the Euclidean Dirac operator in QCD. 
This made it possible to derive exact analytical results for Dirac spectra that have been 
generated by means of lattice QCD Monte-Carlo simulations. Random Matrix Theory 
in this field made it possible to get a much better grip on finite size corrections. 

In Mathematics, Random Matrix Theory was applied to the correlations of the non- 
trivial zeros of the Riemann ^-function [19]. After the impressive numerical work by 
Odlyzko [20] this has been understood in great detail by Berry and Keating (see [21] for 
a review). More recently, Random Matrix Theory has been applied to growth processes 
[22, 23], the longest increasing subsequence of random permutations [24] and the theory 
of random partitions [25, 26]. 

There are several books and reviews on Random Matrix Theory that should be men- 
tioned. The classic text is the book of Mehta [4] which emphasizes the orthogonal poly- 
nomial method. The supersymmetric method is discussed in detail in the book by Efetov 
[6]. In the book by Forester [27] the emphasis is on the relation with solvable models and 
mathematical physics. The historical developments of Random Matrix Theory are dis- 
cussed in the comprehensive review by Guhr, Miiller-Groeling, and Weidenmller. [28]. 
A short historical overview of Random Matrix Theory can be found in the special edition 
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of Journal of Physics A which contains a broad selection of recent papers on Random 
Matrix Theory as well as several reviews [29]. Applications to mesoscopic physics are 
reviewed in [30], applications to QCD in [31] and applications to Quantum Gravity in 
[32]. Among the pedagogical reviews we mention a general review of Random Matrix 
Theory [33], an introduction to the supersymmetric method [34, 35] and applications of 
random matrix theory to quantum gravity [36]. 



2. RANDOIM IVIATRIX THEORY 

In this section we give an elementary introduction to Random Matrix Theory. Some of 
the material in the subsections below is based on the encyclopedia article [33]. 



2.1. The Wigner-Dyson ensembles 

The Wigner-Dyson ensembles are ensembles of hermitian matrices with matrix ele- 
ments distributed according to 

P{H)DH = J^e^"^'^^"DH. (1) 

Here, // is a Hermitian N x N matrix with real, complex, or quaternion real matrix 
elements, and j8 is the Dyson index of the ensemble. Its value is one for real matrix 
elements, two for complex matrix elements and four for quaternion real matrix elements. 
Real quaternions can be expressed in terms of the Pauli spin matrices as 

q = ao + iukGk, (2) 

where the are real variables. The measure DH is the product over independent 
differentials. The normalization constant of the probability distribution is denoted by 
The probability distribution is invariant under the transformations 

H^UHU~\ (3) 

where ?7 is an orthogonal transformation for /3 = 1, a unitary transformation for /3 = 2, 
and a symplectic transformations for j8 = 4. This is the reason why these ensembles 
are known as the Gaussian Orthogonal Ensemble (GOE), the Gaussian Unitary Ensem- 
ble (GUE), and the Gaussian Symplectic Ensemble (GSE), respectively. Since both the 
eigenvalues of H and the Haar measure DH are invariant with respect to (3), the eigen- 
vectors and the eigenvalues are independent with the distribution of the eigenvectors 
given by the invariant measure of the corresponding orthogonal, unitary, or symplectic 
group. Notice that eigenvalues of a quaternion real matrix are quaternion scalars and are 
doubly degenerate. 

Exercise. Show that the symplectic ensemble is invariant under symplectic similarity 
transformations. A symplectic transformation is a unitary transformation that leaves the 
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matrix Z invariant where, in quaternion notation, Z is the direct product of i<72 and the 
identity. 

Exercise. Show that the measure, dH, is invariant with respect to the similarity transfor- 
mations (3). Prove this for each value of j8 . 

One way of deriving the probability distribution (1) is to look for a distribution that 
minimizes the information entropy under the constraints that the average, variance and 
normalization of the probability distribution are fixed. The proof goes as follows [37]. 
The sum of the information entropy and the constraints is given by 

S = j DH[-P{H) logP{H) + aP{H) + bP{H)TTH + cP{H)TtH^] . (4) 

Minimizing this with respect to P{H) we obtain 

- 1 - logP{H) +a + bTvH + cTvH^ = 0, (5) 

resulting in a Gaussian distribution for P{H). The maximum entropy approach has also 
been successfully applied to the joint probability distribution of ^-matrix elements [38] 

and the transport properties of quantum dots (see [39, 30]). 

The joint probability distribution of the Wigner-Dyson ensembles is given by 

pi{?ik})Ud^k=^mm'^u^-''^^'^'d^k^ (6) 

k k 

where the Vandermonde determinant is defined by 

A({A}) = n(^,-Az). (7) 

k>i 

The Vandermonde determinant originates from the Jacobian of the transformation 
from the matrix elements as integration variables to the eigenvalues and eigenvectors as 
integration variables. The distribution of the eigenvectors completely factorizes from the 
probability distribution. As an example, let us calculate the Jacobian for j8 = 2. One way 
to do this is to consider the invariant distance 

TrdHdHK (8) 

By differentiating (3) we find 

U-^dHU ^dUA-A5U + dA, (9) 

where dU = U~^dU. For the invariant distance we then obtain 

TndHdH^ = lTn5UK5UK-lTn5UA^5U + TndA^ 

= lY,8Uu5Uli{Ak-Aif + Y.id^kf- (10) 

k<l k 

We can immediately read off the Jacobian from the Lame coefficients and the result is 
given by the square of the Vandermonde determinant. 

Exercise. Calculate the Jacobian for j8 = 1 and j8 = 4. 
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2.2. The chiral ensembles 



A second class of ensembles are the chiral ensembles with the chiral symmetry of the 
QCD Dirac operator [17, 18]. They are defined as the ensembles of N xN Hermitian 
matrices with block structure 

and probability distribution given by 

P(C)Z)C=^det^/(^^t ) ^-^^■■^'^£>C. (12) 

Again, DC is the product of differentials of the independent parts of the matrix elements 
of C, and Nf is a real parameter (corresponding to the number of quark flavors in QCD). 
The matrix C is a rectangular nx (n + v) matrix. The nonzero eigenvalues of the matrix 
D occur in pairs ±Xk. This can be seen as follows. If 

then 



C \ f a \ >\ { a \ /, 

u • (14) 



ct ; V y " V 

Generically, the matrix/) in (11) has exactly |v| zero eigenvalues. Also generically, the 
QCD Dirac operator corresponding to a field configuration with the topological charge 
V has exactly | v | zero eigenvalues, in accordance with the Atiyah-Singer index theorem. 
For this reason, v is identified as the topological quantum number. The normalization 
constant of the probability distribution is denoted by ^ . Also in this case one can 
distinguish ensembles with real, complex, or quaternion real matrix elements. They are 
denoted by j8 = 1, j8 = 2, and j8 = 4, respectively. The unitary invariance of the chiral 
ensembles is given by 

C^UCV~\ (15) 

where U and V are orthogonal, unitary, or symplectic matrices, respectively. For this 
reason, the corresponding ensembles are known as the chiral Gaussian Orthogonal 
Ensemble (chGOE), the chiral Gaussian Unitary Ensemble (chGUE), and the chiral 
Gaussian Symplectic Ensemble (chGSE), in this order. 

Exercise. Show that the eigenvalues of the matrix with structure (11) occur in pairs 
±A or are zero. If C is an n x (n + v) block matrix show that there are exactly v zero 
eigenvalues. 

It is always possible to decompose the matrix C in (12) as 

C=UAV-\ (16) 
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where A is a diagonal matrix with Aj^j^ > 0, and U and V are orthogonal matrices 
for j8 = 1, unitary matrices for j8 = 2 and symplectic matrices for j8 = 4. The joint 
probability distribution for the eigenvalues is obtained by transforming to A, U and V 
as new integration variables. The Jacobian can be obtained in the same way as in the 
previous section and is given by 

j^Y[x!'^^-'Y[\{xl-xi)\^ (17) 

k k<l 

resulting in the joint eigenvalue distribution 

P{{l])d{X] = JV\^{{X^])fY[X^e-''^'^ilUXk. (18) 

k 

where a = j8 — l+j8v + 2Nf. As in the case the Wigner-Dyson ensembles the distribu- 
tion of the eigenvectors and the distribution of the eigenvalues factorize. 

Exercise. Find the Jacobian (17) using only symmetry arguments and dimensional argu- 
ments. 



2.3. The superconducting ensembles 

A third class of random matrix ensembles occurs in the description of disordered 
superconductors [40, 41]. Such ensembles with the symmetries of the Bogoliubov-de 
Gennes Hamiltonian have the block structure 

where A is Hermitian and, depending on the underlying symmetries, the matrix B is 
symmetric or anti-symmetric. The probability distribution is given by 

P{H)DH = J/' exp (^^T^H^H^ DH, (20) 

where DH is the Haar measure and is a normalization constant. For a symmetric 
matrix B the matrix elements of H can be either complex (j8 = 2) or real (j8 = 1). One 
can easily verify that in this case H satisfies the relation 

H'' = TaHTa, (21) 

where 

' -I ) • <22) 

One can then show that if is an eigenvector of H with eigenvalue X, then Fa is an 
eigenvector of with eigenvalue —X. Because H and have the same eigenvalues, 
all eigenvalues occur in pairs ±X. 
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For an anti- symmetric matrix B the matrix elements of H can be either complex 
(j8 = 2) or quaternion real (j8 =4). In that case H satisfies the relation 

= -TsHTs, (23) 

where 

From this relation it follows that if is an eigenvector of H with eigenvalue A, then Ts^ 
is an eigenvector of with eigenvalue —A. Also in this case, all eigenvalues occur in 
pairs ±A. 

The joint eigenvalue probability distribution of these so called Oppermann-Altland- 

Zimbauer ensembles has the same general form as the joint probability distribution 
chiral ensembles given in (18). The value of a = 1 for j8 = 1 and j8 = 4 whereas for 
j8 = 2 the value of a is or 2. 

Exercise. Derive the joint eigenvalue distribution for the Oppermann-Altland-Zirnbauer 
ensembles. 



2.4. Anti-Unitary symmetries and the Dyson index 

The value of the Dyson index is determined by the anti-unitary symmetries of the 
system. If there are no anti-unitary symmetries the Hamiltonian is Hermitian and the 
value of j8 = 2. 

An anti-unitary symmetry operator, which can always be written as A = UK with U 
unitary and K the complex conjugation operator, commutes with the Hamiltonian of the 
system 

[H, UK] = 0. (25) 
We can distinguish two possibilities 

{UKf^l or {UKf = -l. (26) 

The argument goes as follows. The symmetry operator = {UK)^ = UU* is unitary, 
and in an irreducible subspace of the unitary symmetries, it is necessarily a multiple of 
the identity, UU* — XI. Because of this relation, U and U* commute so that X is real. 
By unitarity we have | A | = 1 which yields A = ±1. 

The case (UK)^ = 1 corresponds to j8 = 1 and the case (UK)'^ = — 1 to j8 = 4. In the 
first case it is always possible to find a basis in which the Hamiltonian is real. Starting 
with basis vector 0i we construct y/i — -\-UK^\. Then choose <^2 perpendicular to 
and define V^2 = 02 + UK<j)2. Then 

(02 + t/i^<fc,V^l) 
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{{UK)%,UK\i/i) 

(02,v^ir=o. 



(27) 



The next basis vector is found by choosing 03 perpendicular to xj/i Yi, etc. . In this 
basis the Hamiltonian is real 



The best known anti-unitary operator in this class is the time-reversal operator for 
which U is the identity matrix. 

In the case (UK)^ = —1 all eigenvalues of the Hamiltonian are doubly degenerate. 
This can be shown as follows. If is and eigenvector with eigenvalue X^, then it follows 
from the commutation relations that also U K(j)i^ is an eigenvector of the Hamiltonian with 
the same eigenvalue. The important thing is that this eigenvector is perpendicular to ^k. 



One can prove that in this case it is possible to chose a basis for which the Hamiltonian 
matrix can be organized into real quaternions [42]. The eigenvalues of a Hermitian 
quaternion real matrix are quaternion scalars, and the eigenvalues of the original matrix 
are thus doubly degenerate in agreement with (29). The best known example in this class 
is the Kramers degeneracy for time reversal invariant systems of half -integral spin but no 
rotational invariance. Then the time reversal operator is given by (72K with ((72^)^ = — 1 • 
Later in these lectures we will discuss other anti-unitary symmetries that enter in the 
Dirac operator for non-abelian Yang-Mills fields. They determine the Dyson index of 
the chiral ensembles. 



2.5. Classification of the Random Matrix Ensembles 

The Random Matrix Ensembles discussed above can be classified according to the 
Cartan classification of symmetric spaces [43, 44, 45]. A symmetric space is a manifold 
such that every point is a fixed point of an involutive isometry (i.e. Xi — > — jc,). The 
Riemann curvature tensor is covariantly constant in a symmetric space. A symmetric 
space is best characterized via the notion of symmetric pair. A symmetric pair {G,H) 
is defined as a pair of a connected Lie group G and a closed subgroup H such that an 
involutive analytic automorphism a of G exists with H e Ha, where Ha is the set of 
fixed points of a. Then, with some additional conditions (see the book by Helgason [46] 
for more details) the coset G/H h a symmetric space. 

As an example, consider the group U{p + q) and define Tpq as the diagonal matrix 
with the first p diagonal matrix elements equal to 1 and the remaining q diagonal matrix 



Hki = {Wk^H\ifi) 



= {UKxi/k,UKHxi/iT 
= {xifk^HUKYiT 



(28) 



(29) 
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elements equal to —1. Consider the transformation 



G{g) = rpqgTpq, (30) 

which is an involution because obviously a{o{g)) — g. One can also easily show that 
it is an analytic automorphism. The set of fixed points of a is the group U{p) xU (q). 
Therefore, U {p + q) /U (p) x U (q) is a Riemannian symmetric space. 

A symmetric space can be of the compact, noncompact or the Euclidean type with 
positive, negative or zero curvature, respectively. Each of the random matrix ensembles 
discussed in this section is tangent to one of the large classes of symmetric spaces. The 
complete classification is given in table I [44] where the corresponding symmetric space 
of the compact type is also given. 

TABLE 1. Random Matrix Ensemble, Corresponding symmetric 
space and the value of the Dyson index j3. 



RMT 


symmetric space 


group structure 


P 


GOE 


AI 


U{N)/0{N) 


1 


GUE 


A 


U{N) 


2 


GSE 


All 


U{2N)/Sp{N) 


4 


chGOE 


BDI 


SO{p + q)/SO{p)xSO{q) 


1 


chGUE 


AIII 


U{p + q)/U{p)xU{q) 


2 


chGSE 


CII 


Sp{p + q)/Sp{p)xSpiq) 


4 


AZ-CI 


CI 


Sp{N)/UiN) 


1 


AZ-D 


D 


SO{N) 


2 


AZ-C 


C 


Sp{N) 


2 


AZ-Din 


Din 


S0{2N)/U{N) 


4 



For example, for special unitary matrices we can write U = l+iH -\ . Therefore, 

the Hermitian matrices are tangent to the space A (after dividing out aU{l) factor). As 
another example, the generators U{p + q)/U{p)x U{q) (class All) are given by matrices 
with the structure (11). More examples can be found in the recent review by Caselle and 
Magnea [47]. 



3. SPECTRA OF COMPLEX SYSTEIMS 

In this section we discuss quantum spectra and several frequently used statistics for the 
statistical analysis of eigenvalues. 



3.1. General Remarks 

Typically, when we talk about spectra we think about quantum mechanics. Among 
the best known examples are the spectrum of a harmonic oscillator with eigenvalues 

£„ = (n + ^)^(0, (31) 
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and the spectrum of an angular momentum operator with eigenvalues of given by 

n^j{j+\). (32) 

However, it is the exception that we are able to write down analytical formulas for the 
complete set of eigenvalues of an operator or even a single eigenvalue. One of the 
questions we wish to address is whether we are able to find any generic features in 
spectra even if we do not have any detailed analytical knowledge of the eigenvalues. 

Let me first give examples of different systems for which spectra play an important 
role. One of the most important examples is the Schrodinger equation 

(-V2 + V(^))V^ = AvA, (33) 

of course with boundary conditions on the wave functions. One particularly simple 
class of potentials are billiards where V{x) is zero inside a connected domain and 
infinite outside this domain. In particular, billiards in the two-dimensional plane have 
been studied extensively in the literature. Although it is trivial to solve the Schrodinger 
equation, the complication arises because of the boundary conditions which can only be 
satisfied for a discrete set of eigenvalues. 

A second example is the Helmholtz equation for the electric and magnetic field inside 
a resonance cavity 

(v2 + yt^)^ = 0, (v2 + yt^)5 = 0. (34) 

In this case one can distinguish different types of boundary conditions, transverse elec- 
tromagnetic, transverse magnetic and transverse electric. One case of special interest 
are transverse magnetic modes in a two-dimensional cavity. For we find the wave 
equation 

(v2+/:2)£^ = 0, (35) 

with the boundary condition that vanishes on the boundary. This is exactly the 
Schrodinger equation for a two-dimensional billiard. This equivalence has been ex- 
ploited for the experimental analysis of quantum billiards [48]. 

A third example is that of sound waves in materials. In that case the wave equation is 
given by 

{W^ + k^)u = Q, (36) 

where u can be decomposed in the longitudinal and the transverse displacement. This 
equation is more complicated than either the Maxwell or Schrodinger equation because 
the two types of modes are coupled through reflections at the boundary. Also this system 
has been a valuable tool for the experimental analysis of spectra of complex systems 
[49]. 

A fourth example is the Dirac operator in QCD. In this case the wave equation is given 

by 

y^{d^ + iA^)\if=a\\f. (37) 
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Here, the 7^ are the Euclidean gamma matrices and the are SU (Nc) valued gauge 
fields. We will discuss this example in much more detail later in these lectures. 

As a fifth example I mention the zeros of the Riemann ^ function. This function is 
defined as 

C{s) = Y^^. (38) 

Empirically, it is know that all zeros are on the line Re* = 5 (also known as the 
Riemann conjecture). An intriguing question is what these zeros have in common with 
the eigenvalues of a complex dynamical system (for answers we refer to [21]). 

All differential operators we have discussed up to now are Hermitian operators. They 
can be diagonalized by a unitary transformation and have real eigenvalues, i.e. 

H = UAU-\ (39) 

Later we will see applications of non-hermitian differential operators. They can be 
diagonalized by a similarity transformation with eigenvalues scattered in the complex 
plane. 

In many cases spectra can be observed experimentally. The best-known examples are 
atomic and nuclear spectra. However, a great deal of spectroscopy has been performed 
on molecular spectra as well. As a more recent example I mention quantum dots in 
which electrons are enclosed inside a heterostructure and collide elastically against 
the boundary. Other examples are the spectra of small metallic particles, which are 
determined by the boundary conditions on the surface. In this case the specific heat 
can be related to the spectral properties of the particles. For references and an extensive 
discussion of these and other examples we refer the [28]. 



3.2. Statistical analysis of spectra 

Although it is the exception that quantum spectra can be obtained analytically, it is 
often possible to obtain a large number of different eigenvalues either experimentally 
or by numerical calculations. This cries out for a reduction of information, and it is 
natural the perform a statistical analysis of the spectra. The average spectral density can 
be written as 

p{X) = (^5{X-h)) (40) 
k 

with normalization given by (for a finite Hilbert space) 

jdXp{X)^N. (41) 

The average can be either over an ensemble of different systems (ensemble average) or 
over different parts of a spectrum (spectral average). Of course it is also possible that the 
average is a mixture of both types of averages. 
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One frequently employs the integrated spectral density 

rE 

N{E)^ / p{X)dX. (42) 

J —oo 

Since N{E) jumps by one at the position of each eigenvalue, it is also know as the 
staircase function. 

Typically, the spectral density varies little over the scale of the average level spacing. 
Therefore the average spectral density can be obtained by locally averaging over an 
interval much larger than the average level spacing. Since (p (A)) depends on the specific 
properties of the system one would like to eliminate this dependence. This is achieved 
by the so-called unfolding procedure. The unfolded spectrum is given by 

K"'= [^\pi^))dX. (43) 
Jo 

One can easily verify that the average spacing of the unfolded sequence is equal to unity. 

The statistical analysis of the unfolded eigenvalues can be performed by a variety of 
statistics. The simplest statistic is the so called nearest neighbor spacing distribution. 
It is denoted by P{S) and is just a histogram of neighboring unfolded levels. A second 
class of statistics is obtained from counting the number of levels in q intervals along 
the spectrum that contains n levels on average. If the number of levels in consecutive 
intervals is given by we can define the moments 

Mp{n) = -j^4. (44) 

Of course, we have that Mi (n) — > n for a large sample. The number variance is defined 
by 

Y?{n)=M2{n)-n^. (45) 

Similarly, one can define higher order cumulants. 
Another frequently used statistic is the A3 statistic which is defined as 



A3(L) ^^j\L'-2L'r + r')Z\)dr. 



(46) 



It has the property that quadratic functions are projected to zero by the kernel of this 
integral operator. Generally, this statistic is much smoother than the number variance 
and for that reason it has been used widely in the literature. 

Exercise. Calculate the number variance and the A3 statistic for a picket fence spectrum, 
i.e. a spectrum of equally spaced eigenvalues. 



3.3. Statistics for uncorrelated eigenvalues 

Let us calculate P{S) and E2(n) for independently distributed eigenvalues with av- 
erage level spacing equal to unity. Then P(S)dS is the probability that there are no 
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eigenvalues inside the interval [0,5] and one eigenvalue in [5,5 + ^^5']. Dividing the first 
interval in n equal pieces we find 

P{S)dS=(l-^ DS^e-^dS for n^oo. (47) 



The number variance can be expressed as 



Y?in) = {rp{X)dXrpa')dX')-n^, 
Jo JO 

= r rdXdX'(T5{X-Xk)5{X'-Xi))-n^, 
Jo Jo Tl 



(48) 



where the Xk are the unfolded eigenvalues. The average factorizes for different eigenval- 
ues. For a finite number of levels we find 

l2(„) = r r dXdX'[5iX-X') + ^''^~^^ '£(^{X-Xk)){d{X' -Xi))]-n\ 
Jo Jo k I 

(49) 

The first term in the brackets is the diagonal contribution with k = I, and in the second 
term we have corrected the absence of the diagonal terms in the sum by a factor 
N{N — l)/N^. After integration over X and X' we obtain 

L2(n) = n-^, (50) 

where A^^ is the total number of eigenvalues. 

The A3 -statistic is obtained by a simple integration resulting in 

A3{L)^^ + 0{l/N). (51) 

However, if one calculates P{S) and £2(0) for spectra of interacting systems one 
finds different results. In particular, one notices a suppression of small spacings and 
a strongly reduced number variance for large n. These properties are characteristic for 
spectra of many complex systems and are known as level repulsion and spectral rigidity, 
respectively. 



3.4. Correlation functions 

The number variance can be expressed in terms of the two-point correlation function 
of the unfolded eigenvalues. The two-point correlation function of the not necessarily 
unfolded eigenvalues is defined by 

P2{X,X') = (p(A)p(AO) - {piX)){p{X')). (52) 
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The nearest neighbor spacing distribution cannot be expressed in terms a two-point 
correlator only. The reason is that it measures the probability that none of the other 
eigenvalues are inside the interval [0, S] . 

The correlation function P2 (A , A') includes a term in which the eigenvalues are equal, 
and can thus be decomposed as 

P2(A, A') = 5(A - A')(p(A)) +/?2(A,A'). (53) 

The two-point correlation function satisfies the sum rule 

JdXp2{X,X')=0, (54) 

where the integral is over the complete spectrum. 

In the literature one also frequently uses the quantity 72 (A, A') for the two-point 
correlation function of the unfolded eigenvalues [4] 

F2(A,A') = -i?2(A,A'), (55) 

where the minus sign is conventional. If a smoothened average spectral density exists, it 
is natural to expect that Y2 is translational invariant, i.e., 

y2(A,A') = 72(A-A'). (56) 
In that case, the number variance can be expressed as 

I?{L) =L- [\l - r)Y2{r)dr. (57) 
^0 



4. THE SUPERSYMMETRIC METHOD IN RANDOM MATRIX 

THEORY 

The supersymmetric has been introduced in Random Matrix Theory and the theory of 
disordered systems to calculate quenched averages according to [7, 6, 50] 

(58) 

z'=z 

The determinant can be expressed as a fermionic integral and the inverse determinant 
as a bosonic integral. For z' = z the generating function has an exact supersymmetry. 
In this lecture we will discuss in detail the supersymmetric method. We start with the 
introduction of Grassmann variables and Grassmann integration. Our main objects are 
graded vectors and graded matrices that consist of both bosonic and fermionic variables. 
As application we will discuss the one-point function and the two-point function of the 
Gaussian Unitary Ensemble. For more details we refer to the book by Efetov [6] and [9] 
where the class j8 = 1 is worked out in detail. 



1 \_:. / det(£> + z) \ 
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4.1. Definitions 

Grassmann variables are anti-commuting variables 

to,Z/} = 0, (59) 

which are necessarily nilpotent, i.e. Xk = 0- can introduce two types of conjugation, 
namely, conjugation of the first kind with the property that 

X**=X, (60) 
and conjugation of the second kind with the property 

X** = -X- (61) 

Below it will become clear that it is often advantageous to use conjugation of the second 
kind. 

The integration over Grassmann variables is defined by 



Jdx = 0, JxdX = ^, 



(62) 



and thus more resembles a differentiation than an integration. As an example, consider 
the integral 



J dxdx*e = J dxdx*{i -XaX*)=a- 



(63) 



Notice that also the differentials dx are anti-commuting elements. 
Next let us make a change of integration variables, t, = aX- Then 



\ = j^d^= j axd{ax), (64) 



so that we have the relation 



d{ax) = -dx- (65) 
a 

For a general transformation given by = ^kiXi we have that 

d^v d^n = -T-rdXx ■ Xn- (66) 
detA 

Exercise. Prove this relation. 

An important class of integrals are Gaussian Grassmannian integrals. Let us consider 

1 = jdxv dXndxl ■ ■ ■ dX*ne-X'^^'Xi . (67) 
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Because of the nil-potency and the definition of the Grassmann integration we obtain 

1 = jdxv --dZudxt ■ ■■dx:^{-XkAkiXir- (68) 

By renaming the integration variables one easily finds that n\ of the terms contributing 
to the integral are the same. We thus find 

/= / ^^;:i---t/;:„^if;:r---^/;:n*(-i)"£;:iAi;,(i);:*(i)---;:„A„;,(„);^;(„). (69) 

where the sum is over all permutations of 1 By rearranging the differential it easily 
follows that the result of Grassmann integrations is exactly the sign of the permutation, 
and 

/ = detA. (70) 
This should be compared to a bosonic integral 

Exercise. Prove this result. 

An essential difference between fermionic and bosonic integrals is that fermionic inte- 
grals are always convergent, whereas bosonic integrals only exist for specific integration 
contours. 



4.2. Graded vectors and Graded Matrices 



An (m|n) graded vector is defined by 



Sm 



(72) 



The matrices acting on these vectors have the structure 



M 



a o 
p b 



(73) 



where a and b are matrices with commuting entries and p and o are matrices with anti- 
commuting entries. The block structure is such that a is an m x m matrix and & is an 
nxn matrix. 
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The trace of a graded matrix is defined by 



TrgM = Tra-Tr&. (74) 

A graded trace of a product of graded matrices is invariant under cyclic permutations. 
The graded determinant is defined by the relation 



detgM = expTrglogM = expTrglog 



a \ r-. f a 
b i^^ + U ip 



(75) 



The logarithm can be expanded in a power series and only the even powers contribute to 
the graded trace. Resumming the power series again into a logarithm, one finds 

detgM= ^ (76) 
deto 

Exercise. Show that detg(A5) = detgAdetg5. 
The transpose of the matrix M is defined by 

^'={-lr^)- (77) 

Because of the additional minus sign we have the relation 

{MiM2f =M^M{. (78) 

The Hermitian conjugate is defined in the usual way 

M'f=M*^. (79) 

A super-matrix is Hermitian if M'^ = M. For a scalar product the Hermitian conjugate 
satisfies the usual relation 

(01,M(|)2)-(M"f01,02), (80) 

where we have used the conjugation of the second kind. 
The eigenvalues of a supermatrix are defined by 

M(^k = hh- (81) 
For two different eigenvectors of a super-Hermitian matrix M we have 

= (M0^0'). (82) 

This can be rewritten as 

(4-A/)((^fc,0/)=O, (83) 

The Supersymmetric Method in Random Matrix Theory and Applications to QCD Febraary 1, 2008 19 



so that the eigenvectors corresponding to different eigenvalues are orthogonal. As a 
consequence, a super-Hermitian matrix can be diagonalized by a super-unitary matrix, 
where a super-unitary matrix is defined in the usual way 

U^U^l. (84) 

The eigenvalues of a graded matrix with block structure (104) are given by the zeros 
and the poles of the secular equation 

de.g(e-X)= ^^,(^1^^ "^ =0. (85) 

The poles of the secular equation can be obtained from the equation 

As an illustration, it is suggested to the reader to calculate the eigenvalues of a diagonal 
matrix. 

The ordinary parts of the eigenvalues are obtained by putting o and p equal to zero 
and are thus given by the eigenvalues of a and b. The eigenvalues with ordinary part 
equal to one of the eigenvalues of a are given by the solutions of eq. (85), and the 
eigenvalues with ordinary part equal to one of the eigenvalues of b are given by the 
solutions of eq. (86). Notice that eq. (85) is singular for A equal to one of the eigenvalues 
ofb. 

Exercise. Show that the eigenvalues of a (1 1 1) super-matrix 



are given by 

X\ = a + 



_ap_ 
a — b^ 



A2 = b+-^. (88) 
a — b 



Next, let us calculate the following Gaussian graded integral 

^^j^L^rl^'^^'^*''^''''^^^ ^^^^ 

where, for reasons of convergence, the matrix M is a super-Hermitian matrix of the form 
(73). In the block notation introduced earlier in this chapter, the exponent can be written 
as 

rkMkii^i = S^aS + S'^ax + x'pS + x'bx 

= {S'^ + x''pci~^HS + a-^ox)+x\b-pa''^G)x. (90) 
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After shifting the bosonic integration variables, the integral factorizes into a bosonic and 
fermionic piece which we already know how to do. The result is 



j_ det{b-pa ^a) _ detb 



deta det{a- ab~^p) detgM' 



(91) 



The second identity follows by writing the determinant as an exponent of the trace 

of a logarithm, and perform a cyclic permutation of the graded trace. The minus sign 
resulting in the inverse determinant comes from the anti-commutation of the Grassmann 
variables. 

It is now clear why the graded trace was introduced as in (74): A graded Gaussian 
integral is given by the inverse graded determinant. 

Exercise. Calculate the superintegrals for the case that the Grassmannian blocks of M 
are equal to zero. 

Finally, a work of caution about notation. In the physics literature "graded" and 
"super" are used as synonyms. For example, graded vectors and supervectors, graded 
matrices and supermatrices, etc. are synonyms. The same is true for abbreviations. For 
example, Str = Trg, Sdet = detg. We will use both notations interchangeably. In quantum 
field theory, supersymmetry is used for a Grassmannian extension of the Poincare group, 
and is different from the use of supersymmetry in these lectures. 

5. INTEGRATION THEOREMS 

If an integrand is invariant under super-unitary transformations, there exist a number of 
powerful integration theorems. In character they are comparable to a complex contour 
integration. 

5.1. The Parisi-Sourlas reduction 

Let us first consider the simplest case involving the graded vectors 



^ , / = (a*,0*). (92) 



d 

het F{p,p^) be an invariant function, i.e., 

F{Up,U'^p^) = F{p,p^), (93) 

for an arbitrary super-unitary transformation U. In this case the integration theorem 
states that 

^ / dada*dede*F(p,p'^) ^2iF(0,0), (94) 
2n J 
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which is also known as Parisi-Sourlas dimensional reduction [5 1] . To prove this theorem, 
we expand F in powers of the Grassmann variables (see also Appendix L of [9]) 

=Foo + Foir + Fio0 + Fiir0. (95) 

If we choose 

then the r.h.s. is only invariant if Fq\ = Fiq = 0. By performing an infinitesimal transfor- 
mation 

we find that 

Fn = -da*Foo. (98) 
a 

Exercise. Show that Fqq is a function of aa*. 

After performing the Grassmann integration, the integral (94) is given by 

-!- / dada*-da*Foo. (99) 
27t J a 

In polar coordinates the integral can be rewritten as 

-2i r . . , 1 



2^ j pdpd<j)-dpFoo. (100) 

For a function that is well-behaved at infinity we then find that the integral is given by 
2jF(0,0). 

This theorem can be extended to supervectors with an arbitrary, but equal number of 
commuting and anticommuting components. We leave it up to the reader to formulate 
this theorem. 



5.2. The Wegner Theorem 

Next we consider an integration theorem for integrals over invariant functions of a 
supermatrix. This theorem was formulated and proved in lectures by Wegner [52] but 
appeared in published form six years later by different authors [53]. In this case the 
invariant function satisfies the relation 

F{Q)^F{UQU-'), (101) 
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where [/ is a super-unitary matrix. The general form of an invariant function is thus 
given by 

FiQ) = Y,akTrg(f. (102) 

k 

Therefore, automatically, F is also invariant under general super-linear transformations. 
The theorem states that 

JdQF{Q)=cF{X1), (103) 

where A is an arbitrary constant Notice that this constant is not determined for an 
invariant function. The constant c depends on the size of the matrix Q and the choice of 
the integration contour. In the example to be discussed below its value is equal to /. 

We first prove this theorem for the simplest case of a (1 1 1) matrix. In this case, any 
function of the graded matrix 



can be expanded as 



If we choose 



Q={ p I ] (104) 



F{Q) = Foo + oFoi+pFio + opFn. (105) 



, e">' , 



Then, 



Invariance of F{Q) then requires that Fio = Fqi = in the expansion (105). Next 
consider the invariance under 

U=(l ? V (108) 



C 1 



Then to first order in (O and ^ 



-\i;{a-b)+p b + cop + CoJ- ^^"^^ 
From the invariance of F we then find 

{(op + ;o){daFoo + dbFoo) + {b- a) (cop + Cc7)Fii = 0. (1 10) 
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This results in 

daFoo + dbFoo 

Fn = 7 • (111) 

a — b 

Example. Consider the invariant function F{Q) — TrgQ^. In components we obtain 
F(Q) = a^ — b^ + lop, which is in agreement with the above general result. 

After performing the Grassmann integrals we find 

/ DQF(Q) = ±l dadb^-^^«^^. (1 12) 

This integral is not well defined if a and b are along the same path in the complex plane. 
For definiteness let us take a along the real axis and b along the imaginary axis. In polar 
coordinates 

fl = p cos , 

b = ipsm(j). (113) 

Then 

daFoo + dbFoo = i^^dp + ^T^d^ )^oo- (1 14) 

The integral over ^ of the second term gives zero provided that the point p = is suitably 
regularized. The integral over p of the first term is the integral of a total derivative. For 
well-behaved functions F we thus find 

jDQF{Q)^iF{0)^iF{n), (115) 

where the second equality follows from the expansion (102). 

An inductive proof for the general case of this theorem was also given by Wegner 
[52]. He considered the case of an arbitrary (m|n) matrix Q. In that case an integral over 
an {m — n)x{m — n) ordinary matrix remains in the final result (for m > n). However, 
for m = n, the theorem is exactly as in (1 15) but with Q equal to an (m|n) super-matrix. 

The proof of the theorem starts from the observation that an arbitrary (m|n) matrix 
can be decomposed into the (1 1 1) super-matrix 



«2=( r?'' n^'"^" )' (116) 



the (1|1) super- vectors 



^«=( rF'" 1' a = 2,---n + m-l, (117) 
and 

^a = (eal,ea«+m), « = 2, • • -n + m - 1. (118) 
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The remaining matrix elements will be denoted by Qr. We can thus write 



FiQ) = FiQ,P,P,QR). (119) 
We now consider the following invariance 

F{UQU-\UP,PU-\Qr)=F{Q,P,P,Qr). (120) 

and define 

KP,P,Qr) = / DQF{Q,P,P,Qr). (121) 

By using the invariance of the measure one easily obtains that 

I{UP,PU~\Qr) ^I{P,P,Qr). (122) 

We thus can apply the Parisi-Sourlas theorem to the P and P integrations. Next we 
can apply Wegner's theorem for (1|1) matrices to the Q integrations. What remains 
is an integral over the smaller super-matrix Qr, thereby proving the general form of the 
theorem. 



5.3. The Efetov-Zirnbauer theorem 

In general, one is interested in super- integrals of a non-supersymmetric integrand. As 
an example we work out the ( 1 1 1 ) case given by 

J DQ{aoo{Q)+aio{Q)o + aoiiQ)p + an{Q)(yp), (123) 

where Q is parametrized as in (104) and aij{Q) are invariant functions. Because of the 
general decomposition of invariant functions, we immediately find that the integrals of 
the second and third terms are zero. However, as was observed by Zimbauer, there is 
another way of deriving this result. 

The proof by Zimbauer [54] is based on the invariance of the measure. If we apply 
the transformation (109) to the Q variables in the integral we obtain 

jDQaxo{Q)op = J DQaio{Q)ic>+(o{b-a)){p + Cia-b)). (124) 

Since this equation is valid for arbitrary O), we find that the second term in (123) 
integrates to zero. In the same way we find that the third term in (104) integrates to 
zero. 

To better appreciate this theorem let us consider a more complicated integral involving 
an invariant function of Q^Q with Q the (1 1 1) supermatrix given in eq. (104). The matrix 
Q^Q is then given by 

, ( a*a + p*p + Y V (125) 

^ ^ y -ao*+b*p b*b-o*o J ^ ' 
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The eigenvalues of this matrix are given by 



, ^ ^ aao*a + b*bp*p+a*b*ap+aba*p* (a*a + b*b)a*G + p*p 

K\ = a a + p p — 1 7 tttt 5 

a*a — b*b {a*a — b*by 

?i2 = ?ii-a*a + b*b-o*o-p*p. (126) 
A general invariant function of Q^Q can therefore be decomposed as 

FiQ^Q)=Foo + Fn(Jp+F,\(7*p*+F2o(y*(7+Fo2P*p+F22(y*(yp*p. (127) 

where the Fij are functions of the commuting variables only. The invariance is given by 

Q^UQV, (128) 

where U and V are super-unitary matrices. 

Because of this decomposition the super-invariant integral over the product of an 
invariant function and an odd number of Grassmann variables is trivially zero. According 
to a theorem due to Zimbauer [54] any super-invariant integral over the product of an 
invariant function and two Grassmann variables also vanishes. The proof proceeds as in 
our previous example and in this case we consider the infinitesimal transformation 



and 



p b 1 [p b [co 1 



a oY (I -C,\( a o\ ( I -a 
p b ^[-co I [ p b [ -p 1 



(129) 



(130) 



If we apply this transformation to the integrand 

dQF{Q^Q)ao*pp\ (131) 



and use the invariance of F{Q''Q) and the invariance of the measure we find that the 
integral over the product of an invariant function and two Grassmann variables vanishes. 
Since this result was also used in Efetov's calculation [7], this theorem is known as the 
Efetov-Zimbauer theorem. 

One might think that the same argument can be applied to the integral of an invariant 
function. Obviously, because of Wegner's theorem this is not true. The reason is that 
the integral (131) diverges in this case. Typically, the invariant measure diverges if two 
eigenvalues coincide. In the next subsection we will analyze this singularity in detail. 



5.4. Efetov-Wegner terms 

Let us study the following (1 1 1) super-invariant integral 

7 = - / dadbdodpe''^'^^^^, (132) 

71 J 
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where 



The additional factor / has been inserted so that the integral converges for the b- 
integration along the real axis. By expanding the Grassmann variables and performing 
the Gaussian integral one immediately finds that 7=1. Next we choose the eigenvalues 
of Q and o and p as new integration variables. The Jacobian of the transformation from 
(Ai,A2) to (a,&) is given by 



1 + ' 



/ = det( (7*)' .^«-fjp )=/. (134) 

{a-ibf ^ {a-ibf 



We thus find that our integral is given by 

dXidX2dodp 



^ r dAidAzdodp ^_^2^^ ^ ^ 
J Tti ' 



(135) 



where the integral over A2 is along the imaginary axis. However, because of the Grass- 
mann integration the result of the integral is zero. Obviously, something went wrong! 
The only place where we could have made an error is at the point a = ib, where the 
variable transformation becomes singular. 

Let us do the integral by excluding an infinitesimal sphere around this point from the 
integration domain. This approach was followed by Haldane and Zimbauer [55] and by 
Gossiaux, Pluhar and Weidenmiiller [56]. In terms of the original variables this changes 
the integral by only an infinitesimal amount. Our integral is thus given by 



/=lim/e (136) 



with 



le^ J dQe-'^'^Q'eiVa^ + b^-e). (137) 
In terms of our new integration variables we have that 



Va^ + b^^^Xf-Xi- (138) 



and our Q function is given by 



e{\/a^ + b^-e) = e{^Xl-Xl-e)-op ^ 5{^Jxl-Xl-e). (139) 

'xi-xi 



Additional terms due to the 5 -function are known as Efetov-Wegner terms. If we use 
polar coordinates for X\ and ^2//, the anomalous contribution to our integral can be 
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wntten as 



If 1 _ 2 

- I rdrd<^dodp-5{r-e)e pa = 1, (140) 
Tl J r 

which is the correct answer for the Gaussian integral (132. In this case the non- 
anomalous piece is zero. 

In general one is interested in non-superinvariant integrals of the following type 

/ = ldQ{MQ)+paMQ)). (141) 

with /o and /i depending only of the eigenvalues of Q. In terms the new integration 
variables the integral is given by 

1 = Mm + ^J dXidX2fx{MM). (142) 

This result can also be interpreted in terms of the integration measure which is now given 
by 



1 1 dUxl-^) 

—dXidX2dadp + ^, dXidX2. (143) 

6. THE SUPERSYMMETRIC METHOD OF RANDOM MATRIX 
THEORY: THE ONE-POINT FUNCTION 

In this chapter we calculate the one-point function of the GUE by means of the super- 
symmetric method. It is evaluated starting from the RMT formulation of the problem. 



6.1. The One-Point Function 

In this subsection we calculate the average resolvent for the Gaussian Unitary En- 
semble by means of the supersymmetric method. This section closely follows the corre- 
sponding section in the article [33]. The starting point is 



(Giz)) = T7 Tr 



1 / 1 \ Id 



N\ z-H I N dJ 
where the generating function is defined by 



Z(y), (144) 

7=0 



det(z-/f) 
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and the integral is over the probability distribution of one of the random matrix ensem- 
bles defined in the first chapter. The determinants can be expressed in terms of Gaussian 
integrals, 



dot{z-H + J) 



dQt{z-H) 
where the measure is defined by 



[dxi/exJ - (z - H)ki(l)i + ixl {z+J- H)kiXi. 



(146) 



(147) 



For convergence the imaginary part of z has to be positive. 
Exercise. Show that thatZ(O) = 1. 

For the GUE (with Dyson index j8 = 2) the probability distribution is given by 

P{H)DH = ^e-^'^'"^"DH. (148) 

The Gaussian integrals over H can be performed trivially, resulting in the generating 
function 



Z{J) = J dxj/^xp -^Trg 



2N^'^\LjXj^j LjXjXj 



-iY,{(l>*z(^j+X*{z+J)Xj) 
j 



(149) 



where the sums over j run from 1 to A^. The quartic terms in and X can be expressed 
as Gaussian integrals by means of a Hubbard-Stratonovitch transformation. This results 
in 



z(y) 



dy/doexp 



N 
'2 



where 



and 



o = 



Trga2-j£ 
j 



z 
z + J 



X] 



Xj 



(150) 



(151) 



(152) 



The variables Obb and Off are commuting (bosonic) variables that range over the real 
axis. Both Obf and Ofb are Grassmann (fermionic) variables. One easily verifies that 
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by completing the square and carrying out the (7-integrations one recovers the original 
quartic interaction. 

The integrals over the and the x variables are now Gaussian and can be performed 
trivially. This results in the (7-model 



Z(y) = J doexp 



N 



Trg(72-Arrrglog((7 + C) 



(153) 



For 7 = we deal with a super-invariant integral, and, according to Wegner's theorem, 
we have Z(0) = 1. By shifting the integration variables according to a ^ a — ^ and 
carrying out the differentiation with respect to J one easily finds that 



{G{z))^{z-i<yFF). 



(154) 



In the large A'^ limit, the expectation value of Off follows from a saddle-point analysis. 
The saddle point equation for gff is given by 



Off + iz= 1/ Off 



(155) 



resulting in the resolvent 



(156) 



Using the relation between the resolvent and the spectral density we find that the average 
spectral density is a semi-circle. 



6.2. Exact calculation of the One-Point Function 



In fact, the integrals in (153) can be evaluated exactly [57] instead of performing a 
saddle point approximation. If we parameterize G as 




(157) 



the partition function can be rewritten as 

Z(y) = 1^1 dadbdadp ^'^^-^^^^'/^^^^ (158) 

By shifting the integration variables and expanding the Grassmann variables, the integral 
can be rewritten as 

Z(y) = -L /rfaJZ7^^a^^j8^(l-^)(l-A^a/3)e^^((«-^)'+(^+'^+^^^^^^ (159) 
In J a" aib 
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After differentiating with respect to J and performing the Grassmann integrations, the 
resolvent is given by 

G{z) = JaJM^-/^)(^ + ^)^-^««-^)'^^*^''^^'). (160) 
The integrals factorize into known integrals. In particular 

fjxiixre-(-'+''^' = ^Hn{z). (161) 

We thus find that the result of the b integral is real. The integral over a is singular. Its 
calculation was discussed in a paper by Guhr [57]. The first observation is that a contains 
an infinitesimal increment ie. By successive partial integrations we obtain 



J a" {n-\)\J a 



{n-\) 

Using the definition of the Hermite polynomials 



(162) 



^1 e-^M'//„_i((z-a)VAV2) (163) 



and the relation 



Im^— = Tibia), (164) 
a — iE 



we find that 



Im/.aVf(-)' = ^(f)"" .-?^«„_,(zVW2). (165, 

To obtain the spectral density we use that 

iN 

p(A) = — (G(A + /e)-G(A-/e)). (166) 

Liz 

Exercise. Use these results to derive an exact result for the spectral density in terms of 
Hermite polynomials. 



7. THE SUPERSYMMETRIC METHOD OF RANDOM MATRIX 
THEORY: THE TWO-POINT FUNCTION 

In this chapter we discuss the two-point correlation function for the GUE using the 

supersymmetric method. As a warm-up exercise we first consider the replica trick with 
fermionic replicas and bosonic replicas. Contrary to the one-point function we will 
obtain the effective partition function by using symmetry arguments. Both the replica 
calculations and the supersymmetric calculation can be found in [58]. Other useful 
sources are the book by Efetov [6], the papers by Zirnbauer [54, 35], the detailed review 
by Zuk [59] and the lecture notes by Fyodorov [34]. 
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7.1. The Replica Trick 

Before embarking on the super symmetric method, we first study the calculation of the 
two-point function by means of the replica trick. In this case we have two possibilities 
for the generating function: a fermionic generating function or a bosonic generating 
function. The fermionic generating function for the two-point function is defined by 

ZniJuJi) = j dHP{H)dQf{x+Ji + ie+H)def{y+J2-ie+H), (167) 

so that 

{G{x + ie)G{y-ie)) = Vim ^ dj,\j^^Qdj^\j^^Q\ogZn{JuJ2) (168) 

and can be represented as a Grassmann integral. 
The bosonic generating function is defined by 

Z_„(7i,72) = J dHP{H)der''{x + Ji+ie+H)der''{y + J2-ie+H) (169) 

and can be represented as a bosonic integral. The two-point function is also given by 
(168). 



7.1.1. The Fermionic Partition Function 
The determinants in the fermionic partition function can be written as 

det"(x+/i + /e+//)det"(3; + /2-«e+^) = J dxdxoxp[-x'^D^yx'^], (170) 
^here Xk = ixh- ■ ■ ,Xn^xh- ■ ■ ^Xn) and D is the block matrix 

y+J2-i8+H j^'^+y -is+H) (1^1> 

For X = y, Ji = J2 and £ = the symmetry of the generating function is U (2n). One 
easily verifies that the resolvent changes sign when e crosses the real axis so that for 
X = 3; the expectation value 

{X'[*x'i) = -{xtx^)- (172) 

This implies that the U (2n) symmetry is broken spontaneously to U{n) x U (n) . In the 
large limit the massive modes can be integrated out and we are left with a theory of 
Goldstone modes. The Goldstone fields are in the coset Q ^ U {2n)/U{n) x U (n) and 
are coupled to the microscopic variables as XkQkiXi- Therefore, they transform as 

Q^UQU^' (173) 
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under a transformation of the microscopic theory as 



Xk^UkiXl with UeU{2n). (174) 

The U (2n) invariance is broken by the mass term. However, invariance can be recovered 
if we transform the mass matrix according to 

M^UMU'K (175) 

The low energy-effective theory should have the same transformation properties. To first 
order in M we can write down only one invariant term, 

TrMQ, (176) 

resulting in the effective partition function 

Zn= [ rfQe^^TrMG ^^^^^ 

jQeU{2n)/U{n)xU{n) 

Here, dQ is the U {2n) invariant measure to be discussed below. We will consider only 
the case that jc and y are near the center of the spectrum. In particular, we take x — (o/2 
and y = — (o/2 so that the mass matrix for Ji —J2 — 0is equal to 

M= with L3 = ) • (178) 

The constant L is determined by the requirement 

lim dj, \j,=j2=o-^, logZ„ = {G{x)), (179) 

which is equal to Kp (0) /N for x = 0. The partition function (177) is the exact low energy 
limit of (167) in the thermodynamic limit for fixed (oN. 
We can parameterize Q as. 



Q^ULiU^ with U 



exp 



A 




(180) 



where A is an arbitrary complex matrix. In a polar decomposition of A — vJ^Ovj with 
vi eU {n) and V2&U{n)/U"{l) and a positive definite diagonal matrix, the matrix U 
is given by 

vj ™* ^=-9.^=-e. (181) 



Therefore, Q can be written as 

v^^ J \ -li^X J \ V2 
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The integration measure in (177) is given by 

Yl[SUn]ki[SU2i]ki, with dU = U-^dU. (183) 

The Jacobian from [dUi2\ki and [5U2i]ki to the variables [Svij^^/ and [5v2]w is given by 
the determinant (for k^l) 

detf -'^f' '^l^' \=Xl-Xl (184) 

For k = I the transformation from [5i2];tA: and [52i]A:/fe to the new variables [5vi]jtfe and A.^ 
so that the Jacobian is given by 

detf -'/,^')^2Xk. (185) 

The integration measure for the new variables is therefore given by 

Yl[5UnUSU2i]ki = Yl^hll(^k-^if- (186) 

Using Uk = A^ — 11^ as new integration variables we finally obtain 

with the thermodynamic limit taken at fixed 

x = N(0Z. (188) 

The partition (187) can be written as a T function. The first step is to expand the 
Vandermonde determinant 

J —t ^ /TTT 



where a and n are permutations of { 1 , • • • , n}, and sg denotes the sign of the permutation. 
Next we use that 



f Yldukufe'''"^ = (a«)«Zi(x), (190) 



-1 

■ k 

which results in 

Znix) = n![det(^&)'+^Zi(;c)]o<ij<„-i. (191) 
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7.7.2. The Bosonic Partition Function 



The treatment of the bosonic partition function is quite similar to the fermionic 
partition function. Now the determinants in the partition function can be written as 

det-"(;c + yi + /e+tf)det~"(y + y2-/£+7/) = ^^nf^n j #i^^02exp[/0;;^D^v0v], 

(192) 

where ^k— ; " " " ; ; 02 ' ' ' ' ' *^*2 ) ^'^^ ^ the block matrix 

The essential difference with the fermionic case is that we have to make sure that the 
integrals are convergent. This is the reason that the sign of the lower right block of 
D has been reversed. For x = y, J\ = J2 and £ = the symmetry of the generating 
function is not U{2n) but rather U{n,n). As in the fermionic case, this symmetry is 
broken spontaneously to U{n) xU (n). In this case, we have that for x = 3; 

(E0f>f)=+(1 02*02)- (194) 

k k 

In the large limit the massive modes can be integrated out by a leading order saddle 
point approximation and we are left with a theory of Goldstone modes. The Goldstone 
fields are in the coset Q E U {n,n) /U {n) x U (n) and transform as 

Q^UQU-' (195) 

under a transformation of the microscopic theory as 

^k^Uu^h (196) 

where U e U{n,n) satisfies 

U%U = Z3. (197) 

In this case the Goldstone modes couple to the microscopic variables as ^^L^Qki^i- The 
U {n, n) invariance of the microscopic theory is broken by the mass term. However, this 
invariance is recovered if we transform the mass matrix by 

M^UMU-\ (198) 

The low energy-effective theory should have the same transformation properties. To first 
order in M we can only write down the term 

TrMQ, (199) 
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resulting in the effective partition function 

Z_„ = / t/e^'^TrMQ (200) 

jQeU{n,n)/U{n)xU{n) 

The mass matrix is defined as in the fermionic case and the constant L is also determined 
by the requirement 

lim --^o'y, |y,=y2=ologZ_„(x = 0) = {G{x = 0)) = 7cp{0)/N. (201) 

n— »0 J\n 

In the second identity we have used that Im(x) > 0. We can parameterize Q as 



Q^UL3U ' with t/ = exp(^j^^ 



(202) 



where A is an arbitrary complex matrix. Notice that there is no factor / in the exponent 

in the definition of U. In a polar decomposition of A = Vj^^0V2 with vi G U{n) and 
V2 e [/(«)/[/"( 1) and jl a positive definite diagonal matrix, the matrix U is given by 



^=('r v'0(m 0(o 'J 



(203) 



Therefore, Q can be written as 



Also in this case the integration measure is given by (186). The Jacobian from [5Ui2]kt 
and [dU2i]ki to the variables [5vi]jt/ and [5vi]jt/ is given by the determinant (for k^l) 

detf"Y' ^i^' ]=^i-^k- (205) 
For ^ = / the new variables are [^vij^ and so that the Jacobian is given by 

The integration measure is therefore also given by (186). 

Using Uk = + jU^ as new integration variables we finally obtain 

poo 

Z-nix)= YldukYliuk-uife'''^^"^ (207) 

•^1 k k<l 

The bosonic partition function can also be written as a T-function. By expanding the 
Vandermonde determinant we can express this generating function as a determinant of 
derivatives 

Z-n{x) = n! [det(a«)'+^Z_i]o<y<n-i . (208) 
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7.2. The Supersymmetric Generating Function 



Next we discuss the supersymmetric generation function for the two-point function 
defined by 



2(71,72) = J dHP{H) 



det(x + 7i + ie + H) dQi{y + J2- ie + H) 
det(x - 7i + ie + H) dQi{y - 72 - ie + H) ' 



(209) 



Here, e is an infinitesimal increment, and P{H) is the Gaussian probabiUty distribution 
for the unitary ensemble. The correlation function of two resolvents is given by 



1 



{G{x + ie)G{y - ie)) = ^ dj, dj, |_,^^oZ(7i ,72) . 
The generating function can be written in terms of Gaussian integrals 

Z(7i,72) = J dHP{H) J d^rJ^r*e''t>tM+'^+")'^^+'Xt{x-Jl+i£+H)Xl 

X {y+J2-ie+H)<j>2+iX2 {y-h-i£+H)xi ^ 

where the vector ^ is defined by 



(210) 



(211) 



/ 01 \ 

\ X2 J 



(212) 



and 



d\l/=f\ ^[^'^^2m#1m#2m 

^l=l 



7t 



The matrix D is defined by 



D = 



/ x+Ji + ie+H 

x-Ji+ie+H 


V 



LM + 







-\y+J2-ie+H] 





(213) 




y-J2-ie+H J 

(214) 



Two widely used matrices that enter in the generating function are the matrix L, which 
breaks the symmetry between the 1 -space and the 2-space, and the matrix k which breaks 
the supersymmetry. They are defined by 



/ 1 
10 

0-10 

V 1 



k = 



/ 1 \ 
0-100 

10 

V -1 / 



k2 

k2 



(215) 
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As in the previous sections we introduce the energy difference 



G) = x-y. (216) 

For (0 = and J\—J2 — the partition function (209) is invariant under the transfor- 
mations 

Xif^Uxif, \i/*^\i/*U'^ (217) 

with 

U^LU=L. (218) 

This group of superunitary transformations is known as ?7(1, 1|2). This symmetry is 
broken explicitly by the mass term ii/^Mij/ and spontaneously by the expectation value 
of 



{xi/*LuXl/i)=NI.T3u. (219) 



Here, Tt, is defined by 



73 



/ 1 \ 
10 
0-10 

V -1 y 



(220) 



This result follows from the same analysis as in the case of bosonic or fermionic replicas. 
The symmetry of the partition function is thus broken spontaneously to?7(l|l)x?7(l|l). 
The Goldstone manifold is given U{l,l /2) /?7(l|l)x?7(l|l) and can be parameterized 
as 

Q = UT3U-\ (221) 

with U G [7(1,1 1 2). The symmetry transformation (217) becomes a symmetry of the 
microscopic partition function if we also transform the mass term by 

M^UMU~\ {111) 

Example. The mass term is given by i//*ZJWi//. Under i// — t/)^ it transforms as 

\if*U^LMU\if=\if*U^LUU-^MU\if. (113) 

Since U^LU = L we recover invariance if M — > [/ MU . 

Because the microscopic fields and the Goldstone fields are coupled as \i/*LQ\f/, the 
Goldstone fields transforms as 

Q^UQU \ (224) 
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To lowest order in M we thus can write down only one nontrivial term with the same 
transformation properties and the microscopic generating function: 



TrgMQ. (225) 

Notice that 2^ = i. 

As before we will study the two-point function in the center of the spectrum. For zero 
sources the mass matrix is given by 

M=|r3. (226) 

We are interested in the thermodynamic limit at fixed CON. To leading order in the 
partition function is given by the integral over the saddle-point manifold, 

2(71,72) = f jQ^i^ZTrg{m+J,iNLTrgk2Qn+J2iNLTrgk2Q22 (227) 

yeec/(i,i|2)/(/(i|i)x(/(i|i) 

with the Goldstone manifold can be parameterized as as in (221). The supersymmetry 
breaking matrix k2 was introduced in (215) and is given by by 

^2 = ( \y (228) 

The subscripts 1 1 and 22 refer to the projection onto the respective blocks of the Q 
matrices. By differentiation with respect to the source terms we thus obtain 



{G{x+ie)G{y-i8)) = / /eTrg(^<2n)Trg(^Q22).^^Trg(r3e)_ 

(229) 



We have seen that the boson-boson block of Q is noncompact. This results in a mass with 
a positive real part when Q is expanded in the generators of U. Because of the supertrace 
we have an additional minus sign for the fermion-fermion term which would result in a 
negative real part. To compensate this minus sign we have to include an extra factor / in 
the generators of U so that the fermion-fermion block becomes compact. The integral 
is therefore over the maximum Riemannian submanifold ofG/(2|2)/G/(l|l) xG/(l|l) 
and the measure is the invariant measure of the unbroken group [44, 60]. 

If X and y would be on the same side of the real axis, the matrix would not appear 
in the partition function, and the symmetry between the 1 -variable and the 2- variables 
would not be spontaneously broken. In that case there is no Goldstone manifold and, 
to leading order in l/N the integrals over the 1 -variables and 2-variables factorize. 
Therefore, they only contribute to the disconnected piece of the correlation function. 



7.3. The Integration Manifold 

In order to identify the integration manifold, we will be guided by the principle of 
convergence. Grassmann integrals are always convergent, and we will ignore them for 
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the moment. Then fermion-fermion blocks and the boson-boson blocks decouple and 
can be treated separately. This corresponds to the case of one fermionic replica or one 
bosonic replica discussed earlier in this chapter. 

In the case of one bosonic replica for x and one for y, the symmetry transformation 
satisfies 

U^a3U = 03 (230) 
with 03 the third Pauli matrix. If we parameterize U as 



we obtain the equations 



a*a — b*b = 1, 
c*c — d*d = — 1, 
ac* = bd*, 
CO* — db* . 

(232) 

Combining these equations we find that \a\ = \d\ and \b\ = \c\.ln our partition function, 
the matrix U only enters in the combination of the coset U1.3U^^. Without loss of 
generality we can omit the phase factors of a and d. Then we obtain b = c*. Our result 
for the parameterization of U is thus given by 

U=( ^"n^ /r^y (233) 

Next we consider the fermion-fermion block which enters in the effective partition 
function as 

e>-^[TrG«^-Tre^^]. (234) 

With some foresight we have chosen the signs of the Hamiltonian the same in the two 
fermionic integrals in the generating function. The symmetry of the fermionic variables 
is therefore C/(2). The Goldstone manifold is still given by Qpp = UcJsU'^. Because U 
is now compact the expansion of Qff to second order in the generators of U results in 
a convergent integral. If we would have started with a noncompact formulation also for 
the fermionic variables the integral over Qpp would not have been convergent because 
of the additional minus sign due to the supertrace. In fact, we would have been forced to 
rotate the integration contour of Qff so that we again would have arrived at a compact 
parameterization of Qff . 

Convergence of the integrals thus forces us to parameterize U according to the com- 
pact subgroup of G/(2) which is C/ (2). Therefore, Qff e C/(2)/C/(l) x C/(l), which can 
be parameterized as 



COS0 «e"^sin0 
je~"^sin0 COS0 



(235) 
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More suggestively this can be written as 



We observe that the parameterization of the fermion-fermion variable is structurally the 
same as the parameterization of the boson-boson variables. 

Next we consider the Grassmann variables. A matrix U G i7(l,l|2) contains 4 
Grassmannian entries. However, in the coset U four of them can be absorbed by the 
?7(1|1) x?7(l|l) subgroup. Our final result for the parameterization of U is then given 
by 

with 

= ( ^ S ) ' ^21 = J* ) , (238) 

and ot, and 5 Grassmann variables. From (218) we find a relation between the 
Grassmann variables by conjugation. However, because a Grassmann variable and its 
conjugate are independent integration variables, the conjugation relation can be ignored 
as has been done in the above parameterization. 



7.4. The Invariant Measure 

Now that we have found a parameterization of the integration manifold we are able 
to calculate the invariant measure. From measure theory we know that the invariant 
measure is given by 

dctg l^^^^ytndtn (239) 

with 5U = U"^dU . This measure is invariant under ?7(1, 1 12) and left invariant under 
[/(1|1) xC/(l|l) [60]. 

Instead of calculating the graded determinant for the general case we hope to convince 
the reader that it is equal to unity by evaluating it for two special cases. The general case 
is left as an exercise. 

Example. We consider the bosonic matrix 

\ a V 1 + « « / 
In order to calculate the invariant measure on the space of U matrices we have to evaluate 

l(a*da — add*) \/l + a*ada — '^^'^^'^ +a da) 



U^^dUi = 



2^ ^ 2Vl+a*fl 

ad a +u'' da) 1 , 

2y/\+a*a 2* 



\/l +a*ada* — " '■"'y' '^"^ i {ada* — a* da) 
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(241) 



The Jacobian for the transformation from the variables {U^ ^dUi)i2 and {U^ ^dUi)2i to 
the variables a and a* is given by the determinant of the matrix 



_*2 



(242) 



which is easily seen to be equal to unity. The invariant measure is thus simply given by 
dado*. 

Example. Next we consider a more complicated example with U2 given by 



U2 = 

For U2^dU2 we obtain 
/ -^{ad5 + 5da) 



( 1 + ^ 


5 



1 + 



i3y 





i8 



r 




1 + 



7)8 



a 





(243) 



1 + f / 



V 







d5 





-i(7^/j8 + j8J7) 






dp 

-^(ydP + Pdr) 








■j{add + dda) J 



(244) 



One easily sees that the Jacobian for the transformation from the 12 and 21 blocks of 
1/2^ dU2 to a, j8, 7 and 5 is equal to unity. The invariant measure is therefore given by 
dad^dydS. 

Exercise. Calculate the Berezinian for the general case and show that its value is equal 
to unity as well . The invariant measure is thus given by dti2dt2\. 

In terms of the tn variables the measure is flat with invariant line element given by 



ds =Trgdti2dt2i. 
This measure is invariant under the transformations 



(245) 



ti2 Ut\2V~ 



til —>■ Vt2\U 



(246) 



7.5. Transformation to Polar Variables. 

Although an explicit calculation of the Grassmann integrations using the parame- 
terization of the previous section is possible, the algebra of such calculation is quite 
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laborious. For that reason we introduce new integration variables defined by the polar 
decomposition 



and 



/Ml \ -1 



. K \ _i 

'21=^1 /Ml r ' 



(247) 



(248) 



with u and v given by 



u = 



1 + ^ a 



V = 



yd 

2 



1 + ^ r 



5 1-4 r 



The matrices u and v satisfy the property that 



U03U = as, 



VO3V = as, 



(249) 



(250) 



with as the third Pauli-matrix. Notice that we have not imposed a conjugation on the 
Grassmann variables. This means that for the Grassmann variables we have relaxed the 
condition [44] U^W^L that defines t/(l, 1 12). 

To calculate the Berezinian of the polar variables we consider the variations 



-1 



dtnv — u ^duii+dji — jiv ^dv, 



V dt2\U 



^dvjl* +dil* - il*u ^du 



(251) 



where /i and /i* are diagonal matrices with diagonal elements /ii and 1112, and /i[ 
and respectively. The block of the Berezinian corresponding to differentiation 
of the boson-fermion blocks of u~^5t\2V and v~^5t2\u with respect to jU and jU* is 
equal to zero. The determinant of the block corresponding to differentiation of the 
boson-boson block with respect to these variables is equal to one. The only nontrivial 
contribution to the Berezinian then arises from the variation of the boson-fermion blocks 
of u~^dt\2 and v~^dt2\u with respect to the off-diagonal elements of u~^du and v~^dv. 
The Berezinian from the t\2 and t2\ variables to the polar variables is thus given by the 
inverse determinant of 







V 



-Ml 






Ml 





-Ml 


^'Ml 




\ 

-«M2 


Ml* 



(252) 



The value of the Berezinian is thus equal to 



B 



1 



(Aii*Mi+M2*M2) 



(253) 



The Supersymmetric Method in Random Matrix Theory and Applications to QCD Febraary 1, 2008 43 



To derive the final result for the two-point function we also need 



Trg^2!2ii = Ai+A2 + 2aj8(Ai-A2), 

Trg^2Q22 = -Ai-A2-275(Ai-A2), (254) 
where Ai and A2 are the diagonal elements of Qii given by 

Ai = l + 2|/ii|2, 

A2 = l-2|jU2p. (255) 
For the exponent in the partition function we obtain 

^iNmZ[X,-X2) ^ (256) 

After performing the angular integrations, the final result for the two-point correlation 
function is given by 

X ( Ai + A2 + 2aj8 (Ai - A2) ) (Ai + A2 + 2/5 (Ai - A2) ) . 

(257) 

A factor of 1/ (2;r)^ in the integration measure has been canceled by the angular inte- 
grations. According to the Efetov-Zimbauer theorem, the terms of second order in the 
Grassmann variables do not contribute to the final result. Because of the extra factor 
Ai — A2 there are no problems with singularities. However, the term of zeroth order in 
the Grassmann variables is the product of a divergent result, from the integration over the 
bosonic variables, and zero from the Grassmann integrations. In this case, a careful treat- 
ment of the singularity will result in a finite contribution known as an Efetov-Wegner 
term. 



7.6. Final Result and Efetov-Wegner Terms 

We first consider the term of zeroth order in the Grassmann variables. In section (5.4) 
we discussed the appearance of Efetov-Wegner terms for a Gaussian integral. Following 
the same procedure for the present case, we will regularize the integrals by excluding a 
small region around the singular point a = b = Q. Using the procedure of ref. [56] we 
replace the integration measure by 

dada*dbdb* dada*dbdb* 6 {Va*a + b*b-p), (258) 

and take the limit p — > at the end of the calculation. We express this regularized 
measure in terms of the polar variables jUi, jU^, ^2 and ju|. Some algebra using the 
explicit polar decomposition of tu and ?2i results in 

a* a + b*b = juf jUi + H2H2 + (HiHi - JU2M2) (aj8 + 7^) + 2/jUijU2j87+ 2ifiiH25a 

3 

+ -alirS{lillii+l4l^2), (259) 
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For \/a*a + b*b we obtain 

1 R 
Va*a + b*b = ^ + ttb [(MiMi - /^*M2)(a/3 + yS) +2ijUiju|/37+2/jUi*jU25a] + -ajSyS, 
ZK L 



(260) 



where 



(261) 



Keeping only the terms of zeroth and maximum order in the Grassmann variables the 
Q function can be expanded as 



R 



R' 



e{Va*a + b*b-p) = e{R- p) + -5{R- p)aPr5 + —dR5{R- p)aPYS. (262) 

The derivative of the delta function can be converted into a delta function by partial 
integration. Since our contribution comes from the singularity we only have to take into 
account derivatives of the singular terms. Near i? — > the integrand containing the 5 — 
functions is of the form 



dR 



^5{R-p) + ^dR5{R-p) 



times an angular factor. Partial integration of the second term results in 



dR 



^5{R-p)-^5{R-p) 



^^5{R-p) = ^^5{R^-p'). 



(263) 



(264) 



Performing the Grassmann integration and including the Efetov-Wegner terms (264) 
after transforming back to the ^1,^2 variables in the correlation function (257 results in 



2^2 /.oo /•! 

{G{x+ie)G{y-ie)) = — J dXi j dX^ 



g(Ai-A2-p^) 
2(Ai-A2) 



(Ai+A2)2 + 4 



,jiVraE(Ai-A2) 



(265) 



The first integral can be performed as follows 



dX\ \ dX2 



1 25(Ai-A2-p2) 



-1 



X\ —X2 



(266) 



where we have taken the limit p — > 0. The final result for the correlation function is 
given by 



yz, poo /■ I 

(G(x + ie)G{y -ie))^Y? + — JA2 j ^ JA2e'^^^(^> = + 2i 



.sm{NG)L)e 



iNcoZ 



{coNy 



(267) 



which is the well-known result for the two-point correlation function of the Gaussian 
Unitary Ensemble. The first term is the disconnected piece of the correlation function. 
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8. mXEGRABILITY AND THE REPLICA TRICK 



Integrability and Random Matrix Theory are closely related. In this section we discuss 
the Toda lattice equation and show that the two-point function of the GUE can be 
obtained from the replica limit of the Toda lattice equation [6 1 ] . A closely related method 
which preceded this work is the replica limit of the Painleve equation. For a discussion 
of this approach we refer to the original literature [62]. 



8.1. The Two-Point Function of the Gaussian Unitary Ensemble 

In this section we derive the two-point function of the Gaussian Unitary Ensemble 
from the Replica limit of the Toda lattice equation. In section 7.1.1 we have shown that 
the partition function with n fermionic flavors can written as [63] 

Z„{x) = C„[det(a«)''+^Zi]o<ij<„-i, (268) 

and n bosonic flavors as [61] 

0</,7<n— 1 ) 

(269) 



where 



li{x) = J^^due''"' and Z_i{x) = 



due""*. (270) 



We will now show that the bosonic and fermionic partition functions are related by a 
recursion relation. The main ingredient of this derivation is the Sylvester identity [64, 27] 
which is valid for determinant of an arbitrary matrix A. It is given by 

det(A)Qy;p^, (271) 

where the Qj are cofactors of the matrix A and the Qj-pq are double cofactors. We apply 
this identity for i = j = n — I and / = j = n, and remind the reader that differentiating 
a determinant is differentiating either with respect to all columns or with respect to all 
rows. Because of the derivative structure of the partition function, only differentiating the 
last column or the last row gives a nonzero result, but this is exactly the desired cofactor 
of the matrix with n increased by 1. One can easily show that the relevant cofactors are 
given by 

Cnn — Zfi— I (x) , 

n—\ ~ ^ix'^n— 1 ) 
Cn-ln-l = {dix)'^Zn-l{x), 

(x). (272) 
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By increasing the index by one we obtain the Toda lattice equation 



4logZ„(^)=n \Zn{x)Y' 

where the factor follows from the choice of the normalization constants. We have 
made this choice because the left hand side is proportional to n^. The replica limit of the 
two-point correlation is given by 

1 2 

^1 = - lim ^ dy. logZ„ ix) = Zi {x)Z_ i ix) 

n-»0 

1 /-oo 

1 Ji 

^.sin;ce" ^^^^^ 
= 2/ , (274) 

X X 

which agrees with the exact analytical result for the two-point function. The fermionic 
partition functions, the bosonic partition functions and the super- symmetric partition 
function form a single integrable hierarchy which are related by the Toda lattice equation 
[61]. A closely related way to derive the two-point function of the GUE is to take the 
replica limit and the corresponding Painleve equation. This has the advantage that the 
bosonic partition function does not directly enter in the derivation. For a discussion of 
this result we refer to [62] which preceded our work [61] on the Toda lattice. 



8.2. Toda Lattice Equation 

The Toda lattice equation describes a one dimensional lattice in which neighboring 
atoms interact via a potential that depends exponentially on the distance. The Hamilto- 
nian of such system is given by 

^=2 I I e-(ft-^^-i). (275) 

This integrable system has infinitely many constants of motion. For example, this can be 
concluded from the existence of a Lax pair meaning that the Hamiltonian equations of 
motion can be written in the form 

^ = [^,^], (276) 
where K and L are two operators. Then 

^I^El = nTr[K"- ^ ^] = nTr[K"- ^[K,L]]=0, {111) 
by using the cyclic invariance of the trace. 
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The Hamilton equations of motion can be written as 

= e'lk-^^^k _ gqk-qk+i . (278) 

If we define 

Zk+i=Zke-''^ (279) 

we obtain the recursion relation 

J , ^ Zk+ iZk-i Zk+iZ k 

2 1 



f logzk+i + ^ logz, = ^-^^^^ - (280) 
dr dt^ zi zi 



We then conclude that 



1 Zk+iZk-i , 

—^logzk^ 2 — + ^0- (281) 

at z^. 

The constant cq can be eliminated by a redefinition of Zk — ZA;exp(co?^/2). 

This equation has been studied extensively in the literature and has appeared in 
different areas of physics. As an example, we mention applications to Chem-Simons 
theory [65] and other applications to Random Matrix Theory [66, 27]. The application 
of the Toda lattice equation discussed in the previous section [61] illustrates the power 
of integrability which makes it possible to obtain exact nonperturbative by means of the 
replica trick. More applications along these lines will be discussed later in these lectures. 
Other applications are parametric correlations and the QCD3 spectral density [67] which 
both have been analyzed by the supersymmetric method [68, 69]. 



9. SYMMETRIES OF THE QCD PARTITION FUNCTION 

It is well-known that the QCD action is greatly constrained by gauge symmetry, Poincare 
invariance and renormalizability. These symmetries determine the structure of the Dirac 
operator. In this section we will discuss the global symmetries of the Euclidean Dirac 
operator. They play an essential role for its spectral properties in the deepest infrared. 
In particular, the chiral symmetry, the flavor symmetry and the anti-unitary symmetry of 
the continuum Dirac operator are discussed. 



9.1. The QCD partition function 

The QCD partition function in a box of volume ^3 = L? can be expressed in terms of 
the eigenvalues of the QCD Hamiltonian Ek as 

zQCD = Y,e-^^\ (282) 

k 
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where /3 is the inverse temperature. At low temperatures (/3 o°) the partition function 
is dominated by the lightest states of the theory, namely the vacuum state, with an energy 
density of Eq/V^, and massless excitations thereof. The partition function Zqcd can be 
rewritten as a Euclidean functional integral over the nonabelian gauge fields Api, 

ZQ^^{M)= / JA^]^det(D + m/)e~^^^/^', (283) 
/=i 

where Sym is the Yang-Mills action given by 



Sym = Jd^[~F^y^ + i^^F^^F^ 



(284) 



The field strength is given by 

= a^A« - dyA^^ +fabcA''^A'^. (285) 

The fabc are the structure constants of the gauge group SU{Nc). We also denote the 
gauge fields by 

A^i=A^aY^ (286) 

where are the generators of the gauge group. The dual of the field strength is defined 
by 

Fnv = \e^tvapF''^- (287) 

The integral v = J d'^xF^yF^y is a topological invariant, i.e. it does not change under 
continuous transformations of the gauge field. An important class of field configurations 
are instantons. These are topological nontrivial field configurations that minimize the 
classical action which are classified according to their topological charge V. The pa- 
rameter d in (284) is known as the -angle. Experimentally, its value is consistent with 
zero. In (283) the mass matrix is diagonal, M — diag(mi, ■ ■ -^mNf), but below we will 
also consider a general mass matrix. 

The anti-Hermitian Dirac operator in (283) is defined by 

D = Y^{d^+iA^), (288) 
where the are the Euclidean Dirac matrices with anti-commutation relation 

{7^,7v} = 25^v. (289) 
In the chiral representation the /-matrices are given by 
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In this representation the Dirac operator has the structure 

The integration measure can be defined by discretizing space time 

dAl = YldAl{x). (292) 

X 

A particular popular discretization is the lattice discretization where the QCD action 
is discretized on a hyper-cubic lattice with spacing a. The discussion of lattice QCD 
would be a lecture by itself. For the interested reader we recommended several excellent 
textbooks on the subject [70, 71, 72]. 

Alternatively, we may define the measure in terms of a complete set of orthonormal 
functions 0^^. In terms of the expansion 

= L (293) 

n=l 

the measure is defined as 



Yldan. (294) 

n 

A field theory is obtained by taking the continuum limit, i.e. the limit a^O or N ^ oo 
for the integration measures discussed above. This limit only exists if we simultaneously 
adjust the coupling constant, i.e. g g{a) or g{N). If such limit exists the field theory 
is called renormalizable. For QCD it turns out that this adjusted coupling constant 
approaches zero in the continuum limit, a property known as asymptotic freedom. 

We will be mainly interested in the eigenvalues of the Euclidean Dirac operator and 
how they fluctuate for gauge fields distributed according to the QCD action. We will 
show that below a well-defined scale the fluctuations of the Dirac eigenvalues are given 
by the corresponding chiral RMT discussed in the first chapter. It turns out that QCD 
and chiral RMT have the same low energy limit which we will analyze by means of the 
supersymmetric method. 



9.2. Axial Symmetry 

The eigenvalues of the QCD Dirac operator are given by 

D(j>k^ak(l>k. (295) 

The eigenvalues are purely imaginary because the Dirac operator is anti-Hermitian 
(notice that dji = —d^ and aJ^ = A^). 
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The axial symmetry, or the C/a(1) symmetry, can be expressed as the anti- 
commutation relation 

{75,0} =0. (296) 

This implies that all nonzero eigenvalues occur in pairs ±i^k with eigenfunctions given 
by 0^ and Ys^k- If A-jt = the possibility exists that 75 (^^ ~ (/)^, so that = is an 
unpaired eigenvalue. According to the Atiyah-Singer theorem, the total number of such 
zero eigenvalues is a topological invariant, i.e., it does not change under continuous 
transformations of the gauge field configuration. Indeed, this possibility is realized by 
the field of an instanton which is a solution of the classical equations of motion. On the 
other hand, it cannot be excluded that = while ^^d y^^k linearly independent. 
However, this imposes additional constraints on the gauge fields that will be violated by 
infinitesimal deformations. Generically, such situation does not occur. 

In a decomposition according to the total number of topological zero modes, the QCD 
partition function can be written as 

ZQC0{M) = J^e'^^Z^^^(M), (297) 

V 

where 

^(M) = {ll^}Yli^i + m}))y. (298) 
/ k 

Here, (■ • •) v denotes the average over gauge-field configurations with topological charge 
V weighted by the Yang-Mills action. If we introduce right-handed and left-handed 
masses as complex conjugated masses we find that the 6 dependence of the QCD 
partition function is only through the combination me^^/^f . This property can be used to 
obtain the -dependence of the low energy effective partition function. 

Exercise. Show that the dependence of the QCD partition function is only through the 
combination e^^l^f . 



9.3. Flavor Symmetry 

A second important global symmetry is the flavor symmetry. This symmetry can be 
best explained by writing the fermion determinant in the QCD partition function as a 
functional integral over Grassmann variables, 

ndet(Z) + m/) = / Jv/rfv/e^^'^^/=i'^''('°+'«/)'^'. (299) 
In a chiral basis with i//^ = 75 ^Tr and V^l = — 75 V^l, the exponent can be rewritten as 

Nf 

i^¥{D + ruf) \iff = YrDWr + Wl^Wl + WrMrlWl + WlMlrvL 

(300) 
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where we have used that 



V/^D< = rRYsDYsWi = -rRDwL = 0, (301) 

and similarly for the matrix elements between a left-handed and a right-handed state. To 
better illuminate the transformation properties of the partition function we have replaced 
the diagonal mass matrix by Mrl and Mlr. 

Exercise. Construct i/Ar, 1//^, V^/? and V'l from the original basis. 

For mf = 0'we have the symmetry 

Wr y^R, W^WV'^- (302) 

The only condition to be imposed on U and V is that their inverse exists. If the number 

of left-handed modes is equal to the number of right-handed modes we thus have an 
invariance under GlR{Nf) x Glt{Nf), where Gl{Nf) is the group of complex Nf x Nf 
matrices with nonzero determinant. However, if the number of left-handed modes is not 
equal to the number of right-handed modes, the axial-symmetry group is broken to an 
Sl{Nf) subgroup whereas the vector symmetry with U = V remains unbroken. For m — 
the flavor symmetry is thus broken explicitly to Gly{Nf) x SlA{Nf) by instantons or the 
anomaly. A Glv{l) subgroup of Gly{Nf) corresponds to baryon number conservation 
and is usually not considered when flavor symmetries are discussed. 

What is much more important, though, is the spontaneous breaking of the axial flavor 
symmetry. From lattice QCD simulations and phenomenological arguments we know 
that the expectation value 

(V/v^) = {^rXITr) + {^lWl) ~ -(240Mey)3 (303) 

in the vacuum state of QCD instead of the symmetric possibility {^f^lf) = 0. Phenomeno- 
logically, this is known because of the presence of Goldstone modes. The pions are much 
lighter than the a mesons. The spontaneous breaking of the axial symmetry also follows 
from the absence of parity doublets. For example, the pion mass and the 5 (or uq) mass 
are very different {nin = \35MeV and m§ = 980MeV). 

For fermionic quarks there is no need to extend the symmetry group to GlR{Nf) x 
Gli{Nf). In that case we will only consider the usual SUR{Nf) x SUt{Nf) flavor sym- 
metry and it spontaneous breaking to SUv{Nf). In the case of bosonic quarks we will see 
in section 3.2 that it is essential to consider the complex extension of SU (Nf). Notice 
that by extending the symmetry group to complex parameters the number of conserved 
currents remains the same. 

On easily verifies that {ij/y/) is only invariant for U — V. The vacuum state thus 
breaks the chiral symmetry down to Gly{Nf). Notice that only the axial symmetries are 
broken. This is in agreement with the Vafa-Witten theorem [73] which states that vector 
symmetries cannot be broken spontaneously in vector-like theories such as QCD. We 
also observe that the complete axial group is broken. The reasons behind this maximum 
breaking [74] of chiral symmetry are less well understood. 
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9.4. Anti-Unitary Symmetries and the Three-fold Way 



The QCD partition function with three or more colors in the fundamental represen- 
tations has no anti-unitary symmetries. For two colors and fundamental fermions and 
for adjoint fermions, the Dirac operator has an anti-unitary symmetry which will be 
discussed below. The classification of QCD according to anti-unitary symmetries was 
introduced in [18]. 



by 



9.4.1. QCD with Two Colors 
For QCD with two colors in the fundamental representation the gauge field is given 

k ^ 



where the X]^ are the Pauli matrices acting in color space. Because of the pseudo-reality 
oiSVil) we have that 

Al = -T2A^T2. (305) 
Using the explicit representation for the /-matrices one can easily show that 

f^i = 72747^7274- (306) 
For the Dirac operator iD = iJudn + Yn^n we thus have 

T27274(«^)*1^27274 = 'iD. (307) 

This can also be written as 

[i^Cr5T2,Z)]=0, (308) 

where K is the complex conjugation operator and C = ^74 is the charge conjugation 
matrix. Because 

{KCr5T2f = 1 (309) 

we can always find a basis such that the Dirac matrix is real for any A^. The proof 
goes along the same lines as the proof that time reversal symmetry results in real matrix 
elements for the Hamiltonian in quantum mechanics (see section 2.4). 

Exercise. Construct a basis for which D becomes real (without relying on the properties 
of D other than the anti-unitary symmetry (308) and (309)). 
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9.4.2. QCD in the Adjoint Representation 



For QCD with gauge fields in the adjoint representation the Dirac operator is given by 

D = rn^^,^ab+r''rnAa^„ (310) 

where the denote the structure constants of the gauge group. Notice that 

is antisymmetric in its indices so that D is anti-Hermitian. Because of the complex 

conjugation property of the 7-matrices we have that 

[1^7475^,1)] =0. (311) 

One easily verifies that in this case 

(1^7475^)^ = -1 (312) 

so that the eigenvalues of D are doubly degenerate (see section 2.4). This corresponds to 
the case jS/j = 4, so that it is possible to construct a basis in which the matrix elements 
of the Dirac operator can be organized into real quaternions. 



10. LOW ENERGY LIMIT OF QCD 

In this section we derive the chiral Lagrangian that provides an exact description of QCD 
at low energies. In the domain where the kinetic term of this chiral Lagrangian can be 
ignored this Lagrangian constrains the Dirac spectrum to satisfy sum rules for its inverse 
eigenvalues. 



10.1. The chiral Lagrangian 

For light quarks the low energy limit of QCD is well understood. It is given by a 
chiral Lagrangian describing the interactions of the pseudo-scalar mesons. The reason 
is that pions are Goldstone bosons which are the only light degrees of freedom in a 
confining theory such as QCD. To lowest order in the quark masses and the momenta, 
the chiral Lagrangian is completely determined by chiral symmetry and Lorentz in- 
variance. In the case of Nf light quarks with chiral symmetry breaking according to 
SUtiNf) X SUR{Nf) SUv{Nf) the Goldstone fields are given by f/ e SU{Nf). Under 
an SUtiNf) x SUR{Nf) transformation of the quark fields, 

Vl^UlWl, (313) 
the Goldstone fields U transform in the same way as the chiral condensate 

U^UlUU^K (314) 
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The symmetry (313) is broken the mass term. However, the full syrmnetry can be 
restored if we also transform the mass term in (300) as 

Mrl ^ UrMrlU£\ Mlr ^ UiMlrUr\ (315) 

The low energy effective theory should have the same invariance properties. To second 
order in the momenta and first order in the quark mass matrix we can write down the 
following invariant terms: 

Tr{d^Ud^U^), Tr(M^C/), Tr(Mzj?C/-i). (316) 

Since the QCD partition function is invariant under Mr^ ^ Mlr, the effective partition 
function should also have this symmetry. The action of the Goldstone fields is therefore 
given by the so called Weinberg Lagrangian [75, 76] 

=^eff(f/) = ^ Tvid^Ud^U"^) - I Tr{MRLU+MLRV), (317) 

where F is the pion decay constant, and E is the chiral condensate. The second constant 
can be identified as the chiral condensate because it is the mean field value of the mass 
derivative of the partition function. The Goldstone fields can be parametrized as 

C/ = exp(/^n„?7F), (318) 

with the generators of SU(Nf) normalized according to Trt^t^ = 5^. This chiral La- 
grangian has been used extensively for the description of pion-pion scattering ampli- 
tudes [76]. 

To lowest order in the pion fields we find for equal quark masses 

^es{u) = \d^n^d^n^ + ^n«n«. (3i9) 

This results in the pion propagator 1 / {p^ + m^) with pion mass given by the Gellmann- 
Oakes-Renner relation 

ml = ^. (320) 

In the long-wavelength limit the order of magnitude of the different terms contributing 
to the action of the chiral Lagrangian is given by [76] 

S = j d'x^{U)r.L'^-^n%2^^ + L'^{U-zM^ + U%2^^). (321) 

Here, represents the zero momentum modes with no space time dependence, 
whereas the nonzero momentum modes are denoted by by n^2M(^) • decomposition 
has two immediate consequences. First, for 'Lm/F'^ 3> l/V the fluctuations of the pion 
fields are small and it is justified to expand U in powers of IT". Second, for 

— « — (322) 
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the fluctuations of the zero modes dominate the fluctuations of the nonzero modes, and 
only the contribution from the zero modes has to be taken into account for the calculation 
of an observable. In this limit the so called finite volume partition function is given by 
[77] 

Zf(M, d)^ [ ^u^Vm.TrMUe''^'f ^ (323) 



uesu{Nf) 



where the 0-dependence follows from the dependence of the QCD partition function on 
the combination me'^/^f only (see section 9.2). We emphasize that any theory with the 
same pattern of chiral symmetry breaking as QCD can be reduced to the same extreme 
infrared limit. 

The effective partition function at fixed V follows by Fourier inversion 

1 f^^ 

Zv(M) = — / dde-'''^Z(M,e). (324) 
2n Jo 

Combining the integral over SU (Nf) and the integral over U{1) and using that 

detUe'^/^f = e'^, (325) 

we find that 

zf{M)= f det^(C/)e^^'^«TrMt/_ (326) 

JU{Nf) 

We also could have derived this partition function directly from an invariance argu- 
ment. For field configurations with topological charge V the difference of the number of 
right-handed modes and left-handed modes is also equal to v. Under the transformation 
(313) and the transformation of the mass matrix (315), the QCD partition in the sector 
of topological charge V thus transforms as 

ZjcD ^ det^(C/LC/^-i)Z5cD, (327) 

where Ur and Ul are both U (Nf) transformations. To zeroth order in the momentum and 
first order in the quark mass, this covariance property immediately gives the partition 
function (326). 



10.2. Leutwyler-Smilga Sum Rules 

The Leutwyler-Smilga sum-rules [77] are obtained by expanding the partition func- 
tion at fixed V in powers of m before and after averaging over the gauge field configu- 
rations and equating the coefficients. This corresponds to an expansion in powers of m 
of both the QCD partition function and the finite volume partition function (326) in the 
sector of topological charge v. As an example, we consider the coefficients of in the 
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sector with v = 0. We expand the effective partition function to second order in the mass 
and use the group integrals 



DUUf. = 0, 

UeU{n) ^ 

1 



DUUijU^i' = -5ii5jk. (328) 

One easily verifies that the second integral is consistent with the unitarity relations for 
U. We then find 

Zf{m) = 1 + ^m2. (329) 

On the other hand, the expansion of the QCD partition function in powers of m is given 
by 

= 1 (330) 
Equating the coefficients of results in the Leutwyler-Smilga sum-rule [77] 

where the prime indicates that the sum is restricted to nonzero positive eigenvalues. 
Exercise. Show that in the sector of topological charge v the sum rule for Nf massless 

flavors is given by 

1 y2T/2 



By equating higher powers of one can generate an infinite family of sum-rules 
for the eigenvalues of the Dirac operator. However, they are not sufficient to determine 
the Dirac spectrum. The reason is that the mass in the propagator also occurs in the 
fermion-determinant of the QCD partition function. However, as will be shown in the 
next chapter, the Dirac spectrum can be obtained from a chiral Lagrangian corresponding 
to QCD with additional bosonic and fermionic quarks [78]. We conclude that chiral 
symmetry breaking leads to correlations of the inverse eigenvalues which are determined 
by the underlying global symmetries. 



11. THE DIRAC SPECTRUM IN QCD 

In this section we give an interpretation of the Leutwyler-Smilga sum-rules in terms 
of the smallest eigenvalues of the QCD Dirac operator. We show that the smallest 
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eigenvalues of the QCD Dirac operator are related to the chiral condensate by means 
of the Banks-Casher relation. 

The order parameter of the chiral phase transition, ( yy) , is nonzero only below the 
critical temperature or a critical chemical potential. As was shown by Banks and Casher 
[82], (yy) is directly related to the eigenvalue density of the QCD Dirac operator per 
unit four-volume 

L^|(W)|=lim^^M2)), (333) 

where the spectral density of the Dirac operator with eigenvalues {Xj^} is defined by 

p(A) = (£5(A-A,)). (334) 

k 




® 



FIGURE 1. A typical Dirac spectrum. To derive the Banks-Casher relation we integrate the resolvent 
over the rectangular contour in this figure. (Figure taken from [84].) 

Rather than the spectral density we often study the resolvent defined by 

^(^) = E^^' (335) 

which is related to the chiral condensate by 

L = - lim lim G{m)/V. (336) 

The resolvent can be interpreted as the electric field at z of charges at iX^. Using this 
analogy it is clear that the resolvent changes sign if z crosses the imaginary axis. Let 
us look at this in more detail. A typical Dirac spectrum is shown in Fig. 1. The average 
number of eigenvalues in the rectangular contour in this figure is p(A)/. Therefore, if 
we integrate the resolvent along this contour we obtain 

^G{z)^il{G{a + e)-G{a-e))^27j:ip{X)l, (337) 
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where the second identity follows from Cauchy's theorem. Using the symmetry of the 
spectrum we find 



ReG(a + e) = np{?i). 



(338) 



Near the center of the spectrum the imaginary part of the resolvent in negligible which 
immediately results in the Banks-Casher relation (333). 

The order of the limits in (333) is important. First we take the thermodynamic limit, 
next the chiral limit and, finally, the field theory limit. 

An important consequence of the Bank-Casher formula (333) is that the eigenvalues 
near zero virtuality are spaced as 



AX = l/p{0)^7t/LV. 



(339) 



For the average position of the smallest nonzero eigenvalue we obtain the estimate 

Amin=^/SV. (340) 

This should be contrasted with the eigenvalue spectrum of the non-interacting Dirac 
operator. Then the eigenvalues are those of a free Dirac particle in a box with eigenvalue 
spacing equal to AX ~ l/V^^'^ for the eigenvalues near A = 0. Clearly, the presence of 
gauge fields leads to a strong modification of the spectrum near zero virtuality. Strong 
interactions result in the coupling of many degrees of freedom leading to extended states 
and correlated eigenvalues. Because of asymptotic freedom, the spectral density of the 
Dirac operator for large A behaves as VX^. In Fig. 2 we show a plot of a typical average 
spectral density of the QCD Dirac operator for A > 0. The spectral density for negative 
A is obtained by reflection with respect to the j-axis. More discussion of this figure will 
be given in section 12.5. 




FIGURE 2. Schematic picture of the average spectral density of QCD Dirac operator. (Taken from ref . 
[83].) 



Let us study the Leutwyler-Smilga sum rule for equally spaced eigenvalues, i.e. 

nn 



TV 



(341) 
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This results in the sum 



- 1 v2t/2 

E^ = ^. (342) 

n=l ^« " 

We obtain the right parametric dependence of the Leutwyler-Smilga sum rules. The 
density of the smallest eigenvalues is suppressed by the number of massless quarks 
through the fermion determinant. Therefore it is not surprising that the exact result 
decreases with Nf. 

Because the eigenvalues near zero are spaced as ~ 1 /IV it is natural to introduce the 
microscopic variable 

u = Ayr, (343) 
and the microscopic spectral density [17] 

We expect that this limit exists and converges to a universal function which is deter- 
mined by the global symmetries of the QCD Dirac operator. The calculation of this 
universal function from QCD is the main objective of the next chapter. We will calcu- 
late ps{u) both from the simplest theory in this universality class which is chiral Ran- 
dom Matrix Theory (chRMT) and from the partial quenched chiral Lagrangian which 
describes the low energy limit of the QCD partition function. We find that the two re- 
sults coincide below the Thouless energy. Universality can also directly be analyzed 
within the framework of of Random Matrix Theory. This can be done by showing that 
the microscopic spectral density is insensitive to deformations of the probability dis- 
tribution. The interested reader is referred to the extensive literature on this subject 
[85, 86, 87, 88, 89, 90, 91, 92]. 



12. PARTIAL QUENCHING AND THE RESOLVENT OF THE 

QCD DIRAC OPERATOR 

In this chapter we study the spectrum of the QCD Dirac operator by means of the 
resolvent. For simplicity we will only consider the case with three or more colors and 
fundamental fermions. 



12.1. The resolvent in QCD 

In terms of the eigenvalues of the Dirac operator the resolvent is defined as [81] 

1 ^ / 1 



1 f p{X;mi,---,mNf) 



= - / dX ' : (345) 

vJ a+z 
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Here, (■ ■ ■) denotes an average with respect to the distribution of the eigenvalues. Some- 
times the argument z of the resolvent is referred to as the valence quark mass [93]. It 
does not enter in the fermion determinant as a regular quark and is therefore a direct 
probe of the Dirac spectrum. 

The relation (345) can then be inverted to give p(A;mi, . . ..mmj)- As mentioned in 
[78], the spectral density follows from the discontinuity across the imaginary axis, 

Disc|^^,;LG(z) = limG(a + e) -G(a -e) = — J^(5(A + Afe)) = — p(A), 

(346) 

where we have suppressed the dependence on the sea-quark masses. This relation fol- 
lows by writing G{z) as a sum over eigenvalues. 



12.2. Generating function 

In this section we introduce a generating function for the resolvent. Similarly to the 
supersymmetric method discussed before this achieved by introducing the ratio of two 
determinants in the partition function. This corresponds to a Euclidean partition function 
of the form 

zPl(z,z',m/) = jdA [][det(D + m/) e''^^ , (347) 

which we will call the partially quenched or pq-QCD partition function. For z = z' this 
partition function coincides with the original QCD partition function. However, it is also 
the generator of the resolvent (see (345)). In the sector of topological charge v we find 

1 d 



G{z;mi,---,mNf) ^ 



logZPi(z,z',m/). (348) 



In addition to the regular quarks, the partition function (347) has additional bosonic 
and fermionic ghost quarks. Our aim is to find the chiral Lagrangian corresponding to 
(347). If we are successful, we have succeeded in deriving a partition function for the 
extreme infrared limit of the spectrum of the QCD Dirac operator. The derivation of this 
partition function and its domain of validity will be discussed in the next sections. We 
will first discuss the subtleties that arise for the axial symmetry with bosonic quarks. 



12.3. Bosonic Quarks 

For bosonic quarks the Goldstone bosons cannot be parameterized by a unitary matrix. 
The reason is that symmetry transformations have to be consistent with the convergence 
of the bosonic integrals. Let us consider the case of one bosonic flavor. Then 

, _i / m id \ 1 /■ ,2 . ,2 X 



^! \ f m id \ f h 
01 J [id'' m M 02 



(349) 
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so that the exponent is purely imaginary for m = and convergent for Re(m) > 0. 

The most general flavor symmetry transformation of the action in (349) is G/(2) which 
can be parameterized as 

U = e^V with and VV^^l. (350) 

For [/ to be a symmetry transformation for m = we require that 

so that H has to be a multiple of as and V has to be a multiple of the identity. The V 
transformations in (350) are not broken by the mass term and therefore represent the 
vector symmetry. Only the symmetry transformation exp(5C73) is broken by the mass 
term so that the axial transformations are parameterized by 

f/ = ( e-' ) "^^^ (352) 
For Nf bosonic flavors the axial transformations are parameterized by 

U=(^^'^ J^H^ with = H, (353) 

which is the coset Gl{Nf) /U{Nf). 

Also in the case of bosonic quarks the axial symmetry is spontaneously broken: 

(01* 0i) + (0102) = {-^fRXifL) + {WlWr)- (354) 

The reason is that both expectation values are equal to the inverse trace of the Dirac 
operator. 



12.4. The infrared limit of QCD 

Ignoring for the moment convergence questions the global flavor symmetry of (347) 
is given by 

G//f(A^/+l|l) xG/l(A^/+1|1). (355) 

We already have seen that the axial symmetry for bosonic quarks is not U (Nf) but rather 
Gl{Nf)/U{Nf). Although the axial flavor symmetry group of the fermionic quarks is 
not a priori determined by convergence requirements we will see in this section that 
supersymmetry necessarily imposes that this symmetry group is compact, i.e. equal to 

U{Nf). 

Under a GlR{Nf + 1 1 1) x Gl]^{Nf + 1 1 1) transformation the quarks fields transforms 

as 

WR URXj/R, Xj/L UlWl 

Vr^VrUr\ Vl^ VlU[ ^ . (356) 
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The subscripts refer to the right-handed (R) or left-handed (L) chirality of the quarks. 
The latter transformations transform the right-handed and the left-handed fermion fields 
in the same way. 

For M = and V = this is a symmetry of the QCD action. At nonzero mass this 
symmetry can be restored if we also transform the mass term according to 

Mlr^UiMlrUrK (357) 

In the sector of topological charge V the partially quenched partition function transforms 
as 

Zv {Mrl, Mlr) ^ Sdct" UrU^ % {UrMrlU^ i , UiMlrUr ^ ) (358) 

The low energy effective partition function should have the same covariance properties 
as the full QCD partition function. 

Because the expectation value of both the condensate of the bosonic quarks and 
the condensate of the fermionic quarks is nonzero in the chirally broken phase, the 
flavor symmetry of the partially quenched QCD partition function (347) is broken 
spontaneously according to 

GlR{Nf X GlL{Nf + l\l)^Glv{Nf+l\l). (359) 

Therefore the low energy limit of the partially quenched QCD partition function is a 
theory of weakly interacting Goldstone bosons parameterized by 

GlR{Nf + 1 1 1) X GhiNf + 1 \\)/Glv{Nf + 1 1 1) = G/A(iV/ + 1 1 1). (360) 

For a confining theory such as QCD the only low lying modes are the Goldstone modes 
associated with the spontaneous breaking of chiral symmetry. However, the Goldstone 
manifold is not the full group GlA{Nf+ 1|1). As we already discussed for the case of 
bosonic quarks the ^7^(1) transformation in which the bosonic fields with index 1 and 2 
are transformed according to an opposite phase factor is not a symmetry of the partition 
function. What is a symmetry of the partition function is the axial transformation 

<l>2 ^ <fc^^">2. (361) 

Mathematically, this symmetry group isG/(l)/f/(l). Had we restricted ourselves to the 
unitary subgroup U {Nf + 1 1 1 ) of G/ {Nf + 1 1 1 ) from the start, we would have missed this 
class of symmetry transformations. 

Let us now consider the low energy limit of the partially quenched QCD partition 
function in the sector of topological charge V. Under a GlR{Nf+ 1|1) x GlL{Nf+ 1|1) 
transformation the Goldstone modes transform as 

Q^UrQU[\ (362) 
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The low energy effective partition function should have the same transformation prop- 
erties a the partially quenched partition function To lowest order in the mass and the 
momentum we can write down the following invariants 

Sixd^Q-^d^Q and Six{MRLQ) +Six{MLRQ-^). (363) 

If we factorize the Goldstone fields into the zero momentum modes Qq and the nonzero 
momentum modes Q{x) as 

Q = QqQ{x), (364) 

the low energy effective partition function that transforms in the same way as the 
partially quenched QCD partition function (358) is given by 

Zv(M)= / rf(2Sdet^((2o)e"-^^'"^*^^\ (365) 

JQeGl{Nf+\\\) 

where 

^'^{Q) = ^Stid^Q-'d^Q + hstr{MRi^Q) + ^LStr{MLRQ-'). (366) 

We already have seen that the boson-boson block of Gl {Nf -|-l|l)isG/(l)/i7(l).Ifwe 
parameterize the field Q as 

Q = e^kTk^k/F^ (367) 

with the generators of G{Nf + 1 1 1), to second order in the Goldstone fields the mass 
term is given by 




(368) 



Let us take all (ghost-)quark masses positive. Because of the supertrace there is a rela- 
tive minus sign between boson-boson Goldstone modes and fermion-fermion Goldstone 
modes. The boson-boson Goldstone modes are noncompact and require that the overall 
minus sign of the mass term is negative. In order to avoid tachyonic fermion-fermion 
Goldstone modes, we have to compensate the minus sign of the supertrace. This can 
be done by choosing the parameters multiplying the fermion-fermion generators purely 
imaginary. This corresponds to a compact parametrization of the fermion-fermion Gold- 
stone manifold. This integration manifold is the maximum Riemannian submanifold 
[44, 78, 79] of G/ (AT/ + 1 1 1 ) and will be denoted as G/ (A^/ + 1 1 1 ) . 



12.5. Domains in (Partially-Quenched) Chiral Perturbation Theory 

In chiral perturbation theory the different domains of validity where analyzed in detail 
by Gasser and Leutwyler [80]. A similar analysis applies to partially quenched chiral 
perturbation theory [81]. The idea is as follows. The Q field can be decomposed as [80] 

(2 = (2oe'vW. (369) 
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where Qo is a constant (zero-momentum) field. The kinetic term for the \j/ fields is given 
by 



^^^rix)^^^r{x)^L-W{x). (370) 

We observe that the magnitude of the fluctuations of the i/a field are of order 1 /L which 
justifies a perturbative expansion of exp(/i/A(x)). The fluctuations of the zero modes, on 
the other hand, are only limited by the mass term 

^yLStrM((2o + !2o')- (371) 

For quark masses m ^ 1 /VL, the field Uq fluctuates close to the identity and the Uq field 
can be expanded around the identity as well. If m <^ Aqcd we are in the domain of chiral 
perturbation theory. 

The same arguments apply to the partially quenched chiral Lagrangian. However, 

there is an important difference. The mass of the ghost-quarks is an external parameter 
that can take on any value we wish. The mass of the Goldstone modes containing these 
quarks is given by 



2zE 
F2 



M^ = —. (372) 



Therefore, independent of the quark masses there is always a domain where the fluctua- 
tions of the zero momentum modes dominate the fluctuations of the nonzero momentum 
modes. This domain is given by [81] 

z«^_ (373) 

and is known as the ergodic domain of QCD. Sometimes it is also referred to as the 
epsilon regime ^ In order that the non- Goldstone modes do not contribute to the partition 
function we have to require that 

L>-^. (374) 
Aqcd 

In the Dirac spectrum we can distinguish three important scales. First, the scale of the 
smallest eigenvalue, 

Kmn = (375) 



^ The epsilon regime is the regime [80] where m^r ^ 0{£^) and 1/L --^ 0(e) in an expansion in powers 
of e of the chiral Lagrangian. These conditions are more strict than the inequality (373). For example, we 
are still in the ergodic domain if 1/L ~ 0{£^l'^) 
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Second, the valence quark mass corresponding to a valence quark Goldstone boson with 
Compton wavelength equal to the size of the box 



F2 

Wc=^. (376) 

Third, the scale of Aqcd which sets the mass scale of QCD. Based on these scales we can 
distinguish four different domains. In the domain where z is of the order of Amin or less 
we can restrict ourselves to the zero momentum sector of the theory, but we have to take 
into account quantum fluctuations to all orders. Also for Amin ^ z -C mc, we only have 
to include the zero momentum modes but in this case the quantum fluctuations can be 
treated semi-classically, or in field theory language, by a loop expansion. The time scale 
conjugate to rric is of the order of the diffusion time across the length of the box which 
explains the name the ergodic domain. For <IC z <IC ^qcd, chiral perturbation theory 
still applies, but the zero momentum modes no longer dominate the partition function. 
For z ^ AgcD. the masses of the Goldstone modes and the other hadronic excitations 
are of the same order of magnitude. Chiral perturbation theory is no longer applicable 
to the spectrum of the Dirac operator, and one has to take into account the full QCD 
partition function. A schematic drawing of the Dirac spectrum in the different domains 
is given in Fig. 4. 

In the theory of disordered mesoscopic systems the scale below which random matrix 
theory is valid is known as the Thouless energy and is given by [94, 95] 

hD 

Ec = jj, Oil) 

where D is the diffusion constant for the diffusive motion of electrons in a disordered 
sample. The time conjugate to Ec is the time scale over which an electron diffuses across 
the sample. The time scale in mesoscopic physics corresponding to AgcD is the elastic 
scattering time. The domain in between Ec and h/x^, where Tg is the elastic scattering 
time, is known as the diffusive domain. This domain is characterized by diffusive motion 
of electrons in the disordered sample. 

In the derivation of the zero momentum limit of the partially quenched QCD partition 
function we have only relied on symmetries and the spontaneous breaking of symme- 
tries. In particular, this means that any theory with the same symmetry breaking pattern 
and only Goldstone modes as light particles will lead to the same low energy effective 
theory. The simplest theory in this class is chiral Random Matrix Theory introduced in 
chapter 2.2. This theory has the flavor symmetries of QCD and, in the limit of large 
matrices, the non-Goldstone modes decouple from the theory. Therefore, the zero mo- 
mentum properties of QCD can be derived from chiral Random Matrix Theory. 



12.6. Zero momentum limit of the effective partition function 

Since z is a free parameter, independent of the quark masses, we can always find a 
part of the Dirac spectrum inside the domain (373). In this domain the zero momentum 
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limit of the QCD partition function in the sector of topological charge V is given by [78] 



ZUM)^ I dQSdQi^Qe^^^''^^Q+^Q~'\ (378) 

jQeGl(Nf+\\\) 

The number of QCD Dirac eigenvalues that is described by this partition function is of 
the order (see (339)) 

This number increases linearly in A^^ for Nc ^ °° which was studied recently by lattice 
simulations [96]. 

In the next section we will study this partition function in the quenched limit (Nf = 0) 
and show that the resolvent coincides with the result obtained from chiral Random 
Matrix Theory [78, 79]. The Random Matrix result will be derived by means of the 
supersymmetric method in section 12.8. 



12.7. Nonperturbative evaluation of G(z) in the quenched limit 



In this section we evaluate the resolvent of QCD for the simplest case of Nf = and 
V = in the domain z <^ F^/LL^. In this domain the partition function is given by 



V) = I 

Jo 



eeG/(l|l) 



dU exp 



— Str 
2 



z + J 
z 



[Q + Q 



(380) 



where the integration is over the maximum super-Riemannian sub-manifold of G/(l 1 1). 
This manifold is parametrized by [79] 



2 = exp 



a 
j8 



e'1> 
e' 



(381) 



The integration measure is the Haar measure which in terms of this parameterization it 
is given by 



Sdet 



SQ 



ki 



50 dsdadp 



dadfid^ds, 



(382) 



It is straightforward to calculate the Berezinian going from the variables 
{5(2ii, 5(222, 5(2i2, 5(221 } to the variables {50, 5*, 5a, 5j8}. The derivative ma- 
trix is given by 



B = 



52, 



kl 



/ i § 

1 I 

e'-' 

V 



50 5s5a5j8 
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2 ^ 
a 

2 



g-s+i(t) j 



(383) 



Using the definition of the graded determinant one simply finds that Sdet5 = i. Up to a 
constant, the integration measure is thus given by d(j)dsdadp. In general, for Nf ^ 0, 
the Berezinian is more complicated [79]. 
We also need 

= ^^^'-^-l] V (384) 

l}^ ^ ^ y j8(e'^-e-0 cosh5(l-^) ) 

After differentiating with respect to the source term (G(z) = (9/logZ(/) |7=o) this results 
in 

G(z)_ f d(j>dsdadp , ap,^, 

cos0(l+^)— xcoshi(l— ^) 

VZ ~J lit 2 ^ 

(385) 

With the Grassmann integral given by the coefficient of aj8 we obtain 

= /^[cos0+;c(cos0+cosh5)cos0]e^('^°^^-'^°^'^^). 
J 4;r 

All integrals can be expressed in terms of modified Bessel functions. We find [78, 79] 

^ = h{x)KQ{x) + ^{l2{x)KQ{x)+h{x)KQ{x)+2h{x)Kx{x)). (386) 
After using the recursion relation for modified Bessel functions, 

/2(;c)=/o(;c)--/i(:c), (387) 

we arrive at the final result [81, 78, 79] 

G{z) 



x{h{x)Ko{x) +h ix)Ki (x)) (388) 



where x = zVH. 

This calculation can be generalized to arbitrary Nf and arbitrary v. The calculation 
for arbitrary Nf is much more complicated, but with a natural generalization of the 
factorized parameterization, and using some known integrals over the unitary group, 
one arrives at the following expression in terms of modified Bessel functions 

^ ^x{la{x)Ka{x) +Ia+l{x)Ka-l{x)), (389) 

where a = Nf+\y\. This result is in complete agreement with the microscopic spectral 
density derived from chRMT [81]. 

For a = this result is plotted in Fig. 3. We observe that, below some scale, lattice 
QCD data obtained by the Columbia group [97] closely follow this curve. 
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FIGURE 3. The resolvent of quenched QCD. The points represent lattice data obtained by the 
Columbia group, and the theoretical prediction (388) is given by the sohd curve. (Taken from ref. [84].) 



12.8. Random Matrix Calculation of G{z) for the chGUE 

In this section we calculate the one-point function of the chiral Gaussian Unitary 
Ensemble by means of the supersymmetric method [98, 99, 85, 88]. The starting point 
is the partition function 

Z{z,z + J) = Im expf-nL^TrWt], 
^ ' J det(z + D) ^' 

(390) 
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where D is given by the anti-hermitian random matrix 

iW 



In this section we only consider the case of zero topological charge where W is an n x n 
complex matrix. The integral is over the real and imaginary parts of the matrix elements 
of W . For definiteness we take the real part of z to be positive. The resolvent is obtained 
by differentiating the partition function 



G{z)^-T,^ 1 \ _ 1 d\ogZ{z,z + J) 



V \z + D V dJ 



(392) 

7=0 



The determinant in the partition function is written as an integral over Grassmann 
variables x and x*'- 



dct{z + J + D) = Y{dXudXu\[dX2ldX2i 

i=i i=i 



Similarly, the inverse determinant can be written as an integral over the commuting 
variables <j) and 0*: 

dcr\z + D) = Jfld(^ud<t>ufld<hid<hl 



i=l i=l 

'I 



After performing the Gaussian integral over W we obtain 



Z{z,z + J) = lld(l)ud(l)i*dXudXulld(hid(hldX2*idX2i 

•' i=l i=\ 



X exp 



r2 )\Q z )\h \Xl \ Z^J \X2 



X 



exp - ^ {xljXii + 02>h-) {xtiX2j + 01/02^) • (395) 



The terms of fourth order in the fields are decoupled to second order terms by means of 
the Hubbard-Stratonovitch transformation. There are no convergence problems for the 
terms that contain Grassmann variables. In these cases we use the identity 

eXHxxHiln^^ ^cij dpidp2Cxp[-nL\p^ + pi) - {pi-ip2)X*i ■ Xi - {Pi+ip2)X2 ■ X2\ , 

(396) 
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for the term of fourth order in the Grassmann variables and 



(397) 



for the mixed terms. In the above equations, pi and p2 are commuting variables and 
a, a*, j8 and j8* are Grassmann variables. The constants ci and C2 are not necessary 
for the calculation of the partition function. We can always fix the overall normalization 
constant from the requirement that Z(z, z) = 1 . In the case of the the four boson term we 
have to be more careful. This term can be rewritten as 



1 



= exp 



[(01*-01 + 02.0|)2-(0f. 01-02-01) 



*\2i 



(398) 



The uniform convergence of the 0, -integrations necessary to interchange the order of 
the integrations can be achieved if we perform a Hubbard-Stratonovitch transformation 
using the identity [88] 



ni 



/oo poo 
-oo ,/ — oo 



(399) 



where a is real positive and — has a negative real part. Indeed, these conditions are 
fulfilled if we choose 



1 



(0f ■01 + 02-0: 



1 



(0f- 01 -02 -02*). 



(400) 



An alternative way to find correct transformation of the four-boson term is to use the 
Ingham-Siegel integral instead of the Hubbard-Stratonovitch transformation [100, 101]. 
Combining the above results we obtain 

Z(z,z + /) = c J d(l)id(j)id(j)2d(l)2dxidxldX2dX2 d[o]exp[-nZ^{o^ + pf + p^ + a*l5+l5*a 

\ 



X exp- 



01 1 







a 




z + ioe"^ 





xl 


1 a* 







Vzl J 


\ 


j8* 








i8 


Z + P1+P2 ) 



(h 

Xl ' 

\ X2 / 

(401) 



where 



d[G] = odods dpi dp2 da da* dfi d^^ 



(402) 



The Gaussian integrals over the and x variables are equal to the product of inverse 
graded determinants. Writing out the 2 x 2 graded determinants according to (76) we 
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obtain 
Z 



x{z + J + pi-ip2T{z + J + pi+ip2T 
x{z + ioe'+ —) "(z + ioe ' + 



z + J + pi-ip2 z + J + pi+ip2 

(403) 

The Grassmann integrals can be performed trivially by expanding the exponent and 
the pre-exponent. The expression for the partition function can be further simplified 
by shifting crcoshs by iz and pi by —z — J. Using the polar coordinates 

Pi = -z-7 + pcos^/E, 

P2 = psin(^/E, (404) 



0-2p2 



and rescaling <7 by 1 /L we find 

/oo poo poo p2.7t 

odo / ds pdp / d(p 
-oo J-oo Jo Jo 

Because of the shift of the a-variable its integration path lies below the real axis. 

The partition function (405) is valid for any n. In the large n limit the p and the 
o integrations can be performed by a saddle point approximation. The saddle point 
equations are given by 

p2=i, o^ = -l, (406) 

leading to the saddle points 

p = ±l, o = ±i. (407) 

The saddle point p = — 1 and and o = —i are not accessible (remember that o contains 
an infinitesimal increment). By expanding p = 1 + 5p and <7 = 1 + do to second order 
about the saddle point we obtain 

/oo pco poo p2.K 

ddo I ds j ddp / d(p 
-oo J — oo J — oo JQ 

X [2zLcoshs(5(7)2 - 2(z +/)Lcos <p(5p)2 - i(5p^) + ii5of] 

^ ^-2n((5(7)2+(5p)2)+zEcoshi-(z+/)Ecos^»+E2((z+7)2-z2)) ^^^g^ 

The terms that are odd in either dp or 5o integrate to zero and have not been written 
down. The Gaussian integrals over 5p and 5^ are easily calculated. As final result for 
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the partition function we find 



z(z,z+y) = 




X [(z + y)Lcosh5-zLcos<p]e-2«(^'=o^h.-zEcos<p+E2((,+/)2_^)) ^ 

(409) 

The normaUzation constant in (408) has been chosen such that Z(z, z) — I. The integrals 
can be expressed in terms of modified Bessel functions. This results in 

Z(z,z + y) = 2nZ{z + J)Ko{2nLz)h{2nL{z + J)) + 2nZzKi{2nZz)Io{2nL{z + J)).i4lO) 

The normahzed resolvent follows after differentiation with respect to the source term. It 
is given by 




= x{Ko{x)Io{x)+Ki{x)Ii{x)) with jc = 2nZz, (411) 



which we have obtained previously from the low energy effective partition function (see 
eq. (389)). 



13. INTEGRABILITY AND THE REPLICA LIMIT OF LOW 
ENERGY QCD PARTITION FUNCTION 

In Chapter 5 we have shown that the GUE two point correlation function can be obtained 
from the replica limit of the Toda lattice equation. In this section we will show that a 
similar integrable structure exists for the partially quenched QCD partition function. In 
particular, partition functions with bosonic quarks, fermionic quarks and the supersym- 
metric partition function form a single integrable hierarchy related by the Toda lattice 
equation [61, 102, 103]. 



13.1. Partition Function for Nf Fermionic Flavors 

We consider the zero momentum limit of the QCD partition function with n fermionic 
flavors with mass m. The effective partition function for this case was already derived in 
section 10.1 and is given by 

Zf{M) = f det^(C/)^^^'^^TrMf/_ 

Ju{n) 

In a diagonal representation of U this partition function can be rewritten as 

/ YldGkYlle'^'^-e'^'l^e'^'^^^'^'^, (413) 

k<l 
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where x = mVL. By writing the Vandermonde determinant as 

= dtt[e'P%<p<n-l,l<q<n (414) 

k<l 

and expanding the determinant, the angular integrals can be written as modified Bessel 
functions. They can be recombined into a determinant as follows: [104, 105] 

Z^{x) = Cndet[Iv+k~i{x)]o<k,i<n-i. (415) 

Using recursion relations for Bessel functions this determinant can be rewritten as a 
T— function 

^« W = ^^Mixdx)'+'zUx)]o<k,i<n-i. (416) 

with 

Z|'(x)=/v(x). (417) 

The simplest way to prove this is to start from (416). To the columns of the matrix we 
apply the relation 

xd^[xPly+p{x)] = {v + 2p)Iy+p{x) +xP+hy+p+i{x) (418) 

starting with p — and increasing p successively by 1 until all derivatives in the first 
row are gone. The first term in the recursion relation is canceled by the addition of a 
multiple of the previous column. We arrive at 

Z„ = c„det[{xd^)''Iy+i{x)]o<k^i<n-i. (419) 

Next we apply the recursion relation 

xd^[xPly+p{x)] = xP+^Iy+p-i{x) - VxPly+p{x) (420) 

starting with the last row. The second term in the recursion cancels by the addition of a 
multiple of the previous row. We continue until the second row, and after that we start 
again in the same way with the last row and continue until the third row. After repeating 
this procedure until all derivatives are gone we end up with the representation (415) 
which completes the proof. 

Exercise. Work out the proof for n = 2. 



13.2. The Bosonic Partition Function 

We already have seen in section 9.3 that the Goldstone manifold for n bosonic quarks 
is given by Gl{n)/U{n). Using the same invariance arguments as before we obtain the 
low energy effective partition function 

Zr„ = / det^((2)e^^™(e+G-^). (421) 

JQ€Gl{n)lU{n) 

The Supersymmetric Method in Random Matrix Theory and Applications to QCD Febraary 1, 2008 74 



In this case Q can be diagonalized as 

Q = [/diag(e'*)f/-^ (422) 
so that an eigenvalue representation of this partition function is given by 

/ Yl'^SkYle'"'Yl(^''-e''){e-'^-e-'')e''^'^''°'^'K (423) 
k k k<l 

The Vandermonde determinants can again be written as 

Ylie"^ -e")= dct[eP''']o<p<n-i,i<q<n (424) 
k<l 

and a similar expression for —s^. By expanding the two determinants, the integrals 
combine into modified Bessel functions which can be combined into a determinant as 
follows [106] 

Z_n{x) = C-ndct[Ky+k+i{x)]o<k^i<n-i. (425) 

Using recursion relations for Bessel functions this determinant can be rewritten as a 
T— function 

Z_„(x) = -^det[(xa^)^+'z_i(x)]o<fe,/<„-i. (426) 

with 

Z_i{x) = Ky{x). (427) 

This simply follows from the observation that (— l)^^v(^) and Iv{x) satisfy the same 
recursion relations. 



13.3. Replica Limit and Toda Lattice Equation 

Instead of the supersymmetric generating function one could have used the replica 
trick to calculate the resolvent. The fermionic replica trick is defined by 

G(z) = lim-logZ„^(z). (428) 
n— >o n 

and the bosonic replica trick by 

G(z) = lim — logZr„(z). (429) 
n^O —n 

If we take the replica limit of the fermionic or bosonic partition functions directly, 
we will obtain a result that differs from the supersymmetric calculation. We will see 
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below that these problems can be avoided if the take the replica limit of the Toda lattice 
equation. 

We now consider bosonic and fermionic partition functions with all masses equal to z 
which only depend on the combination 

X = zYN. (430) 

Next we use the Sylvester identity [64, 27] which is valid for the determinant of an 
arbitrary matrix A. It is given by 

CijCpq — CiqCpj = det{A)Cij-^pq, (431) 

where the Qj are cofactors of the matrix A and the Qj-^pq are double cof actors. By 
applying this identity to the determinant in (416) for z = j = n — 1 and p = q = n, and 
writing the cofactors as derivatives of the partition function, we easily derive the Toda 
lattice equation [107, 108] 

This equation has also been derived as a consistency condition for QCD partition func- 
tions [109]. 

Exercise. Show that the cofactors that enter in the derivation of (432) can be written as 
derivatives of the partition function. 

In the derivation of (432) we have only used that the fermionic partition function is a 
T-function. Since the bosonic partition function is also a T-function (see (426), it satisfies 
the same Toda lattice equation. 

The two semi-infinite hierarchies are connected by 

\im-{xd,f\ogZ^{x). (433) 
n— »o n 

By extending to Toda lattice hierarchy to include an additional spectator boson, it can be 
shown that [103] 

lim-(;c^^)2logZ^(;c) = \imxd^{xdx + ydy)logZi^^i{x\y) 

= .Ji;^^limj:^;clogZi i(x,3;) 

y—*x 

= xdxxG{x). (434) 

For the resolvent we obtain the identity 

xdxxG{x) = l^Zl (jc)Z^ 1 {x) , (435) 

which explains this factorization property of the resolvent. In the same way it can 
be shown the factorization of the susceptibility in a bosonic and a fermionic partition 
function [102]. 
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Inserting the expressions for Z\ and Z_i we find [61] 

G(jc) = x{Ky {x)Iy {x) + Ky- 1 {x)Iy+i {x) ) + , (436) 

where the integration constant a is fixed by the condition that in the sector of topological 
charge v we have v zero eigenvalues resulting in a contribution of v/xto the resolvent. 
In the calculation with the supersymmetric method we found (436) without the topo- 
logical term because we started from a generating function for the nonzero eigenvalues 
only. This result was first obtained from chiral Random Matrix Theory and has also been 
derived from the replica limit of a Painleve equation [62]. 



14. DIRAC SPECTRUM AT NONZERO CHEMICAL POTENTIAL 

In this Chapter we study the quenched microscopic spectrum of the QCD Dirac operator 
at nonzero chemical potential. Using the replica limit of the Toda lattice equation we 
obtain the exact analytical result [102]. 



14.1. General Remarks 

At nonzero baryon chemical potential the Dirac operator is modified according to 

D^D + ^Yo. (437) 

This Dirac operator does not have any hermiticity properties and its eigenvalues are 
scattered in the complex plane [110, 111, 1 12, 1 13, 1 14]. For small ji we expect that the 
width of the cloud of eigenvalues ~ jU^ [110]. The average spectral density is given by 

p{X) = (^d\X-h)), (438) 

k 

and the average resolvent is defined as usual by 




(439) 



Using that d^*l/z — 7t5^{z) we easily derive 

dz*Giz)\,=^ = ^p{?i). (440) 

In this case it is possible to write down a supersymmetric generating function for the 
resolvent. However, we did not succeed in evaluating this partition function and leave 
it as a challenge to the reader to solve the problem by this method. We will calculate 
the resolvent and the spectral density from the replica limit of the Toda lattice equation 
[102]. 
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The replica limit of the spectral density is given by [115, 116, 117] 



p{z,z*) = lim^dA*^ogZn{z,z*), (441) 
with generating function given by 

Z„{z,z*) = (det«(D + Ai7o + z)det«(-D + Airo + z*)). (442) 

The low energy limit of this generating function is a chiral Lagrangian that is deter- 
mined by its global symmetries and transformation properties. By writing the product of 
the two determinants as [1 18] 

\ 

n ;V/-;/ n r* 

det(D + ;U7o + z)det(-D + jU7o + z*) = det 



id + 11 


z 








id- 11 


z* 


z 


id^ + ju 








z* 


id^- 



y id^ -n J 



(443) 



where we have used the decomposition of the Dirac operator given in (291), we observe 
that the U{2) xU (2) flavor symmetry is broken by the chemical potential term and the 
mass term. Invariance is recovered by transforming the mass term as in the case of zero 
chemical potential and the chemical potential term by a local gauge transformation. 

If we expand the chiral Lagrangian for the spectrum of quenched QCD at nonzero 
chemical potential we obtain the following types of terms 

P^^a, m^«, M^^fl, mLnl (444) 
The zero momentum modes factorize from the partition function if 

^«_, and J2'^J2- (445) 

In this domain the partition function is uniquely determinant by invariance arguments 
as in the case of zero chemical potential. In the addition to the mass term we obtain a 
term from the zero momentum part of the covariant derivative resulting in the partition 
function 

Zn{z,Z*)= I rff/e-^Tr[t/,fi][t;-',iJ] + fTrM(C/+f/-l)^ (44^^ 
JUiln) 



where 



5=1 o' X,)' "=("0 A. I- <^^> 
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14.2. Integration Formula 

By decomposing a U {2n) matrix as 



Ml \ ( v\ \ ( \/\—b^ b 



U2 J \ V2 J \ b -Vl-b^ 



U=( ... )("„)( ^ )( '1 t ), (448) 



^2 



t ' 



with Ml , M2, vi e f/ (n), V2 E U{n) /U^{\) and b a diagonal matrix, the following integra- 
tion formula can be proved [102] 

/ JC/det''C/e2Tr[j»^(f/+f^-']+Ip«pTr[(t/sj/-is)^'] 

Ju{2n) 

= ^^)^„;n-i) ^^t[(xa,)^(H)^Z[(x,y)]o<;^,/<n-l, (449) 

where 

^0 

and Cn is an n-dependent constant. 

Exercise. Show that the Jacobian of of the transformation (448) from U~^dU to 
u~[^dui,U2^du2,vY^dvi,V2^dv2 SiadYikd^k^^ given by 

/= n {^i-^lffli2h), (451) 

\<k<l<n k=l 

where Xk = -y^l — 

Exercise. Show that forn = 1 and all Up—O the integral (jc, _y) is given by 



%/ 

j/v (j:)/v- 1 ( -y) + xiv- 1 (;c)/v ( -y) 



x^ —y'^ 



(452) 



It can be interpreted as the QCD partition function with two different masses [119]. 



14.3. Toda Lattice Equation 

Using the integration formula (449) for p = 1 we find that the zero momentum 
partition function Z^(z,z*) (see eq. (446)) is given by 
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where 

Z\{z,zn= /'ArfAe-2^^V(A2-i)|/^(A^yj:)|2. (454) 
Jo 

By applying the Sylvester identity to the determinant in (453) for / = j = n — I and 
p = q — n and expressing the cofactors as derivatives, we find a recursion relation that 
can be written in the form of a Toda lattice equation 

zd^z d,*logZ^ {z,z) = —{zz ) ^^y^^ ^^^-^2 • (455) 

For the spectral density we find the simple expression (Zq (z,z*) = 1) 

p{zX) = limJ-<9,5,*logZ,^(z,z*) = ^zr(z,z*)zri(z,z*). (456) 
M^o Tin I 

What remains to be done is to calculate the partition function with one bosonic and one 
conjugate bosonic quark which will be completed in the next subsection. 



14.4. The Bosonic Partition Function 



In this subsection we evaluate the low energy limit of the QCD partition function at 
nonzero chemical potential for one bosonic quark and one conjugate bosonic quark. Be- 
cause of convergence requirements, the inverse determinants of nonhermitian operators 
have to regularized [120, 121]. This is achieved by expressing them as the determinant 
of a larger Hermitian operator 



det" 



lim det 



z id + n 
id^ +11 z 

( e 


z* 

\ -id'^ + H 



^ det-1 


/ z* -id + ji 
\ -id^ + ^ z* 





z id + ^ \ 




e 


id^ + ju z 




-id + fl 


e 




z* 


e ) 





(457) 



To better expose the symmetries of this problem we rewrite the inverse determinant in 
the r.h.s. of this equation as 



J d^kd^k ^^p[i 



r2 



( 



\ 



e 

z* 

-id^ + jl 




z 
e 




id + fl 

e 

-z 



\ 

id — jl 

-z* 
e 









02 




03 




V 04 / 



In this case the mass matrices are given by 
Ci = ( * - 1 and 



e z 
z* e 



-z 
e 



(458) 



(459) 
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The two mass matrices are related by 

C2 = -Kil with I={^\1^- (460) 
For jU = and Ci = C2 = the basic symmetry^ of the partition function is G/(2) x G/(2), 

This symmetry can be extended to nonzero mass or chemical potential if we adopt the 
transformation rules 

Ml ^ f/2iUi [/2- \ iU2 ^ f/l iU2[/r ^ , (462) 

where /ii is the chemical potential matrix that is added id and ^2 is the chemical potential 
matrix that is added to —id^. These matrices are introduced for the sake of discussing the 
transformation properties of the partition function (458) and will ultimately be replaced 
by their original values 

^1=^2= _'^)- (463) 

The chiral symmetry is broken spontaneously by the chiral condensate to G/(2) . Because 
of convergence of the bosonic integral, the Goldstone manifold is not G/(2) but rather 
G/(2) /U{2), i.e. the coset of positive definite matrices as in the case of zero chemical 
potential. Under a G/(2) x G/(2) transformation the Goldstone fields transform as 

g^f/2gf/-i. (464) 

The low energy effective partition function should have the same transformation prop- 
erties as the microscopic partition function (458). To second order in and first order in 
the mass matrix we can write down the following invariants 

TrCie, TrC2e"', Tr!2M2e"Vi, TrAiiM2. (465) 

We also have the discrete symmetry that the partition function is invariant under the 
interchange of ^1 and ^2- This symmetry implies that the coefficients of the two mass 
terms in the effective partition function are the same. Using that the integration measure 



This is the symmetry when we disregard convergence. Taking convergence of the bosonic integrals into 
account the vector symmetry is U (2) and the axial symmetry is G/ (2) /f/ (2) . 
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on positive definite Hermitian matrices is given by dQ/die\^Q, we finally arrive at the 
effective partition function^ 

z%(z,z*) = / ^_det-(!2).-^Tr[e,5][e->,B]+f^Tr(Cie+C2e-). 

(466) 

Exercise. Prove that for a parameterization of positive definite 2x2 matrices as Q = 
?7diag(e'*i the invariant measure is given by 

= - e'^)ie-'' - e~'^)dsids2dU, (467) 

det'^Q 

where dU is the invariant measure onU{2)/U{l) xU{l). 

To evaluate the integral (466) we do not use the parameterization of the above exercise 
but rather 



„ t e cosn* e smn* \ /^^o^ 



e'' cosh 5 e'^sinh* 
^ e~'^sinh5 e~''coshs 

where 

re (-00,00), 5 e (-00,00), ^e(_oo,oo), de{0,n) . (469) 

Using that the Jacobian from the measure dQ/det^Q to the measure drdsdtdQ is given 
by 

/ = 4/^ cosh s sinh5, (470) 
we find the partition function {z — x + iy) 

Z\{z,z*) = hmCe [ drdsdtdecoshslsiDhsle^"^ (471) 

e^O J 

Xe'^ (% sinh s cosh t cos 0—4y sinh s sinh t sin +4e cosh r cosh s cosh j) —;U^F-^y( 1+2 sinh'^ s) 



The integral over r is equal to 2^o(2A'^£cosh.vcosh?) with leading singularity given by 
~ — loge. This factor is absorbed in the normalization of the partition function. Then 
the integral over 6 gives a Bessel function. Introducing u = sinhi' as new integration 
variable we obtain [102] 

Zli{z,Z*) = C-ifdt |^^iMe2^%(2VLM(;c2cOSh2f+3;2sinh20l/2)^-M2FV(l+2«2)_ 

(472) 



^ Notice that the Q variables used in this lecture are the transpose the Q variables used in [102]. This 
simplifies the discussion of the transformation properties 
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To perform the integral over u we use the known integral 

£dxx^+'e-^'UPx) = J^^e-P'/'". (473) 

This results in 



(474) 



9 11 /II 

Using that cosh ? = ^ + 2Cosh2? and sinh ? = — ^ + ^cosh2f, the integral over f can be 
rewritten as a modified Bessel function resulting in [102] 

Z\{z^x^iy,z^)^C.,e-^^K, ( ^^^2/2 J • (^75) 



14.5. The Dirac Spectrum at Nonzero Chemical Potential 

The final result for the quenched spectral density is obtained by substituting the 
partition functions Z^(z,z*) and ZX^(z = x + /jjZ*) in expression (456) obtained from 
the replica limit of the Toda lattice equation. We find, 

X f\dle-^^^'^'^'\Iy{lzVZ)\\ (476) 
Jo 

The normalization constant has been chosen such that the jU ^ limit of p {x, y) for large 
y is given by EV/tt (see below). 

In Fig. 4 we show a graph of the spectral density for 3; = and n^F^V = 16 (solid 
curve) plotted as the dimensionless ratio 

^ / N p(^,y) 

p^i^'^y) = -^2^- (477) 

Also shown is the result where the Bessel function Ky is replaced by its asymptotic 
expansion (dotted curve). This result has been obtained from a nonhermitian eigenvalue 
model [122] that is not in the universality class of QCD. An important difference 
between the two results is that the spectral density (476) is quadratic in x for x — > 0, 
whereas the result given by the dotted curve is linear in x for x ^ 0. 

In the limit Re(z)E/ jU^F^ <C 1 the upper limit of the integral in (476) can be extended 
to infinity. Using that 

I Xale-^^-VV-MXzVL)l^ = _pL_,^i^;, (|_^) (478) 
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0.002 



FIGURE 4. The quenched spectral density at nonzero chemical potential in the ergodic domain of 
QCD (full curve). The dotted curve represents the result obtained from a nonhermitian eigenvalue model. 
(Taken from [84].) 



the spectral density can be expressed as 

2 yz^ 

P (-^, y ) = — u^zz*Kv {zz* u)Iv {zz*u) with u = ^^2p2 



(479) 



Therefore, the spectral density becomes a universal function that only depends on a 
single parameter u. This parameter can be rewritten in a more physical way as 



U = 7rpasym(-«,y)- 

For the dimensionless ratio we obtain ^ 

-P^^A^ = 2uzz*Ky{zz*u)Iy{zz*u) 

which is universal combination that depends only on zz*u. 



(480) 



(481) 



This formula arose as a result of a discussion with TUo Wettig. 
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In the thennodynamic limit the Bessel functions can be approximated by their asymp- 
totic Hmit, 

^v(z) = \[fe-\ Iy{z) = -^{e^ + ie-^), Iy{z*) = -y^{e^* - ie'^* ) . (482) 
V 2z y/2nz v27rz* 



This results in 



For y — ;> oo the integral over A can be performed by a saddle point approximation. If 
the saddle point is outside the range [0, 1] we find zero for V ^ oo. This results in the 
spectral density 

and pix^y) = outside this strip. This result is in agreement with the mean field analysis 
of the effective partition function [118]. For the integrated eigenvalue density we find 

/ dxp {x, m) = — (485) 

in agreement with the eigenvalue density at jU = 0. 



15. CONCLUSIONS 

The supersymmetric method is a powerful method to analyze random matrix ensembles. 
This approach makes it possible to formulate the problem in terms of symmetries and 
the spontaneous breaking of these symmetries. The final result for universal correlation 
functions is uniquely given by a partition function of Goldstone modes interacting 
according to an effective Lagrangian with the symmetries of the microscopic partition 
function. Convergence requirements necessary lead to a Goldstone manifold that is 
the product of a compact and a noncompact submanifold. In the final answer for the 
correlation function this structure is present in the form of compact and noncompact 
integrals and is therefore the natural formulation of the problem. 

In these lectures we have given an elementary introduction to Random Matrix Theory 
and the supersymmetric method. We have applied this method to four problems which 
are technically the least demanding: the one-point function of the Gaussian Unitary 
Ensemble (GUE), the two-point function of the Gaussian Unitary Ensemble and the 
one-point function of the chiral Gaussian Unitary Ensemble both as zero and at nonzero 
chemical potential. The one-point function of the GUE is not a universal quantity and 
was only discussed to illustrate the method. The two-point function of the GUE and 
the one-point function of the chGUE are universal quantities. They are relevant for 
eigenvalue correlations of Hamiltonian systems with broken time reversal invariance 
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and no other symmetries and for the distribution of the small eigenvalues of the QCD 
Dirac operator, respectively. 

We have also shown that for j8 = 2 the supersymmetric partition function connects the 
semi-infinite bosonic and the semi-infinite fermionic Toda lattice hierarchy of partition 
functions. This makes it possible to obtain the second derivative of the partition function 
from the replica limit of the Toda lattice equation. The final result is the product of 
a bosonic (noncompact) and a fermionic (compact) integral. In this approach there is 
no need to analyze a supersymmetric non-linear cr-model. The factorized structure of 
the final answer is a direct consequence of the Toda lattice equation and make this the 
minimal approach for the evaluation a supersymmetric partition function. It is still an 
open problem whether integrable structures play a similar role for other values of the 
Dyson index. 

The supersymmetric method is not the only method to derive nonperturbative results 
for Random Matrix Theories. Perhaps the best known method is the orthogonal poly- 
nomial method. This method starts from the joint eigenvalue distribution and exploits 
orthogonality relations to perform the integrals that lead to the correlation functions we 
are interested in. Therefore, this method fails for problems that can not be formulated 
in terms of a joint eigenvalue density. Most notably, the statistical theory of S'-matrix 
fluctuations and the theory of parametric correlations have only been solved by means 
of the supersymmetric method. 

A disadvantage of the supersymmetric method is that it requires a thorough knowl- 
edge of supersymmetry and superanalysis. A straightforward evaluation of the Grass- 
mann integrals leads to a factorial proliferation of terms and does not work in all but the 
simplest cases. To find the right variables to solve the supersymmetric integrals is often 
an arduous task and one wonders why the symmetries of the problem do not provide us 
with a more directed path to the final result. One hint that more progress can be made 
is that in some cases where the supersymmetric method fails because of technical prob- 
lems the problem can be solved exactly by means of the replica limit of the Toda lattice 
equation. We have illustrated this for the quenched spectral density of the QCD Dirac 
operator at nonzero chemical potential. Meanwhile, this problem has also been solved 
by means of the orthogonal polynomial method [123, 124], but the derivation of this 
solution from a supersymmetric nonlinear (7 -model has remained elusive. 
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